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Summary 

The status and some recent developments of 
continuum modeling for large repetitive lattice struc- 
tures are summarized. Discussion focuses on a num- 
ber of aspects including definition of an effective sub- 
stitute continuum, characterization of the continuum 
model, and the different approaches for the genera- 
tion of the properties of the continuum, namely, the 
constitutive matrix, the matrix of mass densities, and 
the matrix of thermal coefficients. Also, a simple ap- 
proach is presented which can be used to generate 
analytic expressions and/or numerical values of the 
continuum properties. 

The proposed approach is applied to some beam- 
like and double-layered platelike lattices, currently 
considered as candidates for large space structures. 
Future directions of research on continuum model- 
ing are identified. These include needed extensions 
and applications of continuum modeling as well as 
computational strategies and modeling techniques. 

1. Introduction 

Lattice structures have been used for many years 
in spanning large areas with few intermediate sup- 
ports. These structures can combine low cost with 
light weight and an esthetically pleasing appearance. 
Also, because of their ease of packaging, transport- 
ing, and assembling in space, lattice structures have 
attracted considerable attention for use in large- 
area space structures such as the Space Station, 
large space mirrors, antennas, multipurpose plat- 
forms, and power systems for supporting space op- 
erations. For the large-area lattice structures con- 
sidered for space applications, a main feature is the 
repetition of a basic pattern or configuration many 
times. 

The state of the art in the analysis, design, and 
construction of lattice structures up to 1976 is re- 
viewed in references 1 and 2. The currently used 
approaches for analyzing large repetitive lattices can 
be grouped into four classes (see table 1), namely, 

1. Direct method 

2. Discrete field methods 

3. Periodic structure approaches 

4. Substitute continuum approaches 

In the first approach (direct method) the struc- 
ture is analyzed as a system of discrete finite 
elements, and the methods of solving structural 
framework problems are applied. It has the obvi- 
ous drawback of being computationally expensive for 
large lattices. This is particularly true when a buck- 
ling, a vibration, or a nonlinear analysis is required. 
Reduction methods could be applied for buckling and 


nonlinear problems to remedy this drawback (refs. 3 
and 4). 

The second approach (discrete field methods) 
takes advantage of the regularity of the structure and 
involves writing the equilibrium and compatibility 
equations at a typical joint of the lattice and either 
solving the resulting difference equations directly or 
using truncated Taylor series expansions to replace 
the difference equations by differential equations (see, 
for example, refs. 5 to 9). This approach works well 
for simple lattice configurations, but becomes quite 
involved for lattices with complex geometry. 

The third group of methods, referred to as 
“periodic structure approaches,” are based on either 

(1) the combined use of finite elements and transfer 
matrix methods, which is efficient only for rotation- 
ally periodic (i.e., cyclically symmetric) structures 
or lattices with simple geometries (refs. 10 to 12), or 

(2) the exact representation of the stiffness of an in- 
dividual member from which the analysis of beam- 
like lattices with simply supported edges can be 
performed (refs. 13 to 15). 

The fourth approach is based on replacing the 
actual lattice structure by a substitute continuum 
model that is equivalent to the original structure 
in some sense, such as the constitutive relations, 
strain energy, and/or kinetic energy (see, for exam- 
ple, refs. 16 to 33). The use of continuum models to 
simulate the behavior of planar lattice beams dates 
back to the previous century (see ref. 34, p. 483). It 
has gained popularity only in recent years and has 
been applied to a variety of other discrete systems 
and phenomena including solid and liquid crystals, 
dislocations and defects, composite materials, and 
biological systems. 

The number of publications on continuum model- 
ing of repetitive lattice structures has been steadily 
increasing. Therefore, there is a need to broaden 
awareness among practicing engineers and research 
workers of the recent developments in various aspects 
of continuum modeling for large lattice structures. 
The present paper is a modest attempt to fill this 
void. Specifically, the objectives of this paper are 

1. To assess the effectiveness of the currently 
used approaches for continuum modeling 

2. To present a simple and rational approach for 
development of continuum models for large 
repetitive lattice structures 

3. To identify the future directions of research 
which have high potential for realizing 

the advantages of continuum modeling. 

The scope of the present study includes thermo- 
elastic stress analysis, buckling, free vibration, and 


geometrically nonlinear problems of large lattice 
structures. Beamlike and platelike repetitive lattices 
with pin and rigid joints are considered. 
Continuum modeling of lattices with flexible joints 
is also discussed. 

2. Advantages of Continuum Modeling 

Before an assessment is made of the different 
approaches for developing continuum models, the fol- 
lowing three advantages of using the continuum mod- 
eling approach for analyzing repetitive lattice struc- 
tures are identified. First, it offers a practical and 
efficient approach for analyzing large lattice struc- 
tures. This is particularly true for beamlike and 
platelike lattices, wherein a dimensionality reduc- 
tion can result in a substantial reduction in the 
number of degrees of freedom. Second, it provides 
a simple means of comparing structural, thermal, 
and dynamic characteristics of lattices with differ- 
ent configurations and assessing the sensitivity of 
their responses to variations in material and geomet- 
ric properties. Third, it provides an effective tool for 
parameter and system identification and for feedback 
control system design of lattice structures. 

3. Definition and Key Elements of a 
Substitute Continuum Model 

A number of definitions have been given for the 
substitute continuum model. Herein an effective 
continuum model is defined to be a continuum with 
the following characteristics: 

1. The same amounts of thermoelastic strain 
and kinetic energies are stored in it as those 
of the original lattice structure when both 
are deformed identically. 

2. The temperature distribution, loading, and 
boundary conditions of the continuum simu- 
late those of the original lattice structure. 

3. For beamlike and platelike lattices the con- 
tinuum models are one-dimensional beams 
and two-dimensional plates, respectively 
(see fig. 1). 

4. Local deformations are accounted for. 

5. Lattices with pin joints are modeled as classi- 
cal continua, and lattices with rigid (and/or 
flexible) joints are generally modeled as 
micropolar continua. 

The last two characteristics are perhaps the most 
important in terms of recent developments and are 
discussed in the next two sections. 

3.1. Local Deformations 

The local deformations of two axially loaded pla- 
nar trusses are shown in figure 2. The first truss has 


double lacing and a single-bay repeating cell. The 
second truss has single lacing and a double-bay re- 
peating cell. The chord members of the first truss 
remain straight as shown on the top sketch. On 
the other hand, the actual deformation of the single- 
laced truss has the zig-zag pattern shown on the top 
right sketch. On the average, however, the chord 
members remain straight. Early continuum models 
averaged these deformations, thereby substantially 
overestimating the axial stiffness. Recent continuum 
models, for lattices with more than one bay in their 
repeating cells, do account for the local deformations 
(see refs. 25, 26, and 35). 

3.2. Ordinary Versus Micropolar Continua 

A contrast between the ordinary and micropolar 
continua is made in figure 3. For an axially loaded 
pin-jointed truss member the transverse motion is 
completely characterized by the joint displacements. 
The member rotation ip is related to the joint dis- 
placements and Wj. Therefore, the appropri- 
ate continuum to use in modeling pin-jointed trusses 
is the ordinary continuum for which the displace- 
ment field completely characterizes the motion of the 
structure. 

On the other hand, for a rigid-jointed member, 
the transverse motion is characterized by both the 
joint displacements w and Wj and the joint rota- 
tions and 8j, which are independent degrees of 
freedom. Therefore, the appropriate continuum to 
use in modeling rigid-jointed flexural members is one 
whose motion is characterized by both a displace- 
ment field and an independent rotation field (referred 
to as “micro- rotation” field). The micropolar contin- 
uum is such a continuum. 

3.3. Characterization of the Substitute 
Continuum Model 

The substitute continuum model is characterized 
by the thermoelastic constitutive relations and den- 
sity parameters, which are determined in terms of 
the geometric and material properties of the origi- 
nal lattice structure. The thermoelastic constitutive 
relations and density parameters of the continuum 
can then be used to determine (1) the thermoelastic 
strain and kinetic energies, (2) the governing differ- 
ential equations, and whenever appropriate (3) the 
equivalent discrete finite element model. 

3.4. Comments on Continuum Models 

The following three comments regarding contin- 
uum models seem to be in order: 
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1. For some lattices the substitute continuum mod- 
els may not have much resemblance to the contin- 
uum theories commonly used in engineering prac- 
tice. Also, for complicated lattices the continuum 
models may be fairly complicated, and therefore, 
not useful for practical applications. 

2. The accuracy of the predictions of the con- tin- 
uum approximation increases with the increase in 
the number of repeating cells (or modules) con- 
stituting the original lattice structure. 

3. The response of the substitute continuum model 
(which simulates that of the original lattice struc- 
ture) can be generated through (1) exact (or an- 
alytic) solution of the governing differential equa- 
tions or (2) application of a discretization tech- 
nique such as the Rayleigh-Ritz technique or the 
finite element method. 

4. Development of Substitute Continuum 
Models for Stress Analysis and Free 
Vibration Problems 

A number of approaches have been proposed for 
developing continuum models and for determining 
the appropriate constitutive relations and density 
parameters. These approaches include 

1. Relating the force or deformation characteristics 
(or both) of a small segment of the lattice to those 
of a small segment of the continuum (refs. 20, 22, 
and 36 to 42) 

2. Using the discrete field method to obtain the 
governing difference equations of the lattice and 
either solving them directly or converting them to 
approximate differential equations (refs. 6 to 8) 

3. Applying homogenization techniques based on us- 
ing multiple-scale asymptotic expansions (refs. 43 
to 46) 

4. Using energy equivalence concepts, the potential 
and kinetic energies of a typical (repeating) cell of 
the lattice are equated to those of the continuum, 
after expanding the nodal displacements of the 
lattice in a Taylor series. 

The latter approach has been applied to a num- 
ber of beamlike and platelike lattices. Computerized 
symbolic manipulation was used to generate analytic 
expressions for the stiffness and density parameters 
of the continuum (see refs. 25, 26, 35, 47, and 48). 
More recently, an equivalent approach was proposed 
for generating the properties of simplified one- and 
two-dimensional continuum models of beamlike and 
platelike lattice trusses with pin joints. This ap- 
proach does not require the use of computerized sym- 
bolic manipulation (see ref. 49). Rather, numerical 
values of the stiffness and mass coefficients can be 
obtained by using a small FORTRAN program on a 


microcomputer. A modified version of this approach 
is now described. 

The three key elements of the foregoing approach 
are 

1. Introduction of kinematic and temperature as- 
sumptions to reduce the dimensionality of the 
continuum 

2. Expansion of each nodal displacement, each strain 
component, and the temperature in a Taylor 
series 

3. Generation of four transformation matrices that 
relate nodal displacements, axial strains, and 
temperatures of individual members of the re- 
peating cell to the displacement, strain, and tem- 
perature parameters of the continuum. 

A schematic representation of the sequence of 
steps involved in the generation of the stiffness coef- 
ficients for simplified one- and two-dimensional con- 
tinuum models is given in figure 4. The proce- 
dure is described subsequently for lattices with pin 
joints, and the appendix lists the computer program, 
which was developed by Sandra L. Whitworth and 
Jeanne M. Peters of the George Washington Univer- 
sity Joint Institute for Advancement of Flight Sci- 
ences, Hampton, Virginia. The procedure consists 
of the three major phases described in the following 
three sections. 


4.1. Phase 1 — Generation of the 
Thermoelastic Stiffnesses of a 
Repeating Cell 

Step 1. A repeating cell (or module) is isolated 
from the lattice grid. The axial strain, temperature, 
and consistent mass matrix of a typical member k of 
the repeating cell are given by e^), J’( fc ) i and [mW], 
respectively. 

Step 2. The axial strain eW is expressed in terms 
of the vector of strain components in the coordinate 
directions through the following matrix equation: 


£<*> = [f#>] { £ }W (1) 


where 


{£}<*> = 


' £ll ' 

£22 

£33 

< > 


712 

713 

• 723 > 


( 2 ) 


m 2 n 2 tm in mn 


(k) 


( 3 ) 
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where £n, £22* and £33 are the axial strains in 
the coordinate directions; 712, 713, and 723 are the 
shearing strains; and (£,m, n)^ are the direction 
cosines of the member. 

For simplicity, in the present study the strain 
state is assumed to be uniform within each repeating 
cell. Variation of the strain state within the repeating 
cell can be accounted for by expanding {£}(*) in a 
Taylor series about the center of the repeating cell. 
The number of terms in the Taylor series expansion 
is equal to the number of independent deformation 
modes of the repeating cell. 

Step S. The stiffness matrix and the thermal load 
vector of the repeating cell are generated as follows: 

[K\= Y2 ( EAL ) {k) [fi®] 1 [fl (fc) ] (4) 

Members 

{P r }= (aEAL) {k '> [iE (fc) ]V fc) ( 5 ) 

Members 

where E, A, L, and a are the elastic modulus, cross- 
sectional area, length, and coefficient of thermal 
expansion of member k; and superscript t denotes 
transposition. 

The thermoelastic stiffnesses of the equivalent 
three-dimensional classical continuum are obtained 
by dividing the right sides of equations (4) and (5) 
by the volume of the repeating cell. Note that for 
members shared by n repeating cells, their cross- 
sectional areas in equations (4) and (5) are divided 
by n. 

4.2. Phase 2 — Generation of Thermoelastic 
Stiffnesses and Density Parameters of the 
Beam and Plate Continua 

Step 4 ■ Reduction in dimensionality is achieved 
by introducing kinematic and temperature assump- 
tions regarding the variation of the displacement 
components, strain components, and temperature in 
the plane of the cross section (for beamlike lattices) 
or in the thickness direction (for platelike lattices). 

The vectors of nodal displacements and strain 
components and the temperature of member k are 
expressed in terms of the corresponding continuum 
parameters by means of transformation matrices as 
follows: 


where {u} c , {e} c , and {T} c are the vectors of dis- 
placement parameters, strain parameters, and tem- 
perature parameters used in describing the beam (or 

plate) continuum; [T^], [T^], and [F ^ ' ] are trans- 
formation matrices. 

Step 5. The stiffness matrix, matrix of density 
parameters, and thermal load vector of the contin- 
uum are given by 

[K\ e = I £ (EAL)W [lf>]‘ [*<*>]* [*<*>] [r<*>] 


Members 


(9) 

w-i e K>] 

[nM] [ri fc) ] 

(10) 


Members 


{Pt}c = ^ E ( aEAL ) {k) [* (fc) ]* |H fc) ] {T}c (11) 
Members 

where Cl is a characteristic geometric property of 
the repeating cell (length of the repeating cell for 
beamlike lattices and planform area of the repeating 
cell for platelike lattices). 

4.3, Phase 3 — Generation of Properties of 
Simplified Beam and Plate Continua 

Step 6. Simplified continuum models are ob- 
tained by partitioning the displacement and strain 
parameters of the continuum as follows: 



<i3) 

where {e c i} are the strain parameters to be retained 
in the simplified model, { e c i } are the strain param- 
eters associated with the neglected stress resultants 
(or internal forces), {e C 3} are the strain parameters 
to be omitted, and {u c i} and {u c 2} are the dis- 
placement parameters associated with the retained 
and neglected mass coefficients, respectively. The 
stiffness matrix, matrix of density parameters, and 
thermal load vector of the continuum are partitioned 
accordingly: 


M W = 

r{ ky 

Me 

(6) 


r*n K12 

K 13 ■ 

{ e }« = 

T { e k) 



[K]c = 

K 22 

K 23 

Me 

(7) 

Symm. 

R 33 . 

y(fc) _ 

rW 
1 T 

{T}c 

(8) 

[•M]c = 

•Mu 

Symm. 

Ml 2 

M22 


(14) 

(15) 
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{Pt)c = I Pt2 1 ( 16 ) 

{ P TZ ) 

The effective thermoelastic coefficients of the sim- 
plified continuum models are obtained by deleting 
the rows and columns associated with {£ c 3} and by 
expressing {e C 2} in terms of {e c \} (using static con- 
densation). The resulting constitutive matrix and 
thermal load vector have the following forms: 

[C} = [K n )~[K 12 }[K 22 }- 1 {K 21 } (17) 

{C T } = {Pri} ~ [Ki2][K22]~ 1 {PT2} (18) 

The density parameters of the simplified contin- 
uum are obtained by neglecting the terms associated 
with {u c 2}, that is, 


M = [Mil] (19) 

Step 7. The thermoelastic strain energy density 
of the continuum can be written in the following 
form: 

Uc^Uo-Ui- U 2 (20) 

where U 0 is the isothermal strain energy density, and 
Ui and U 2 are the contributions to the strain energy 
of the linear and quadratic terms in the temperature 
parameters. The expressions for U 0 and U\ are 

Uo = |{e c i}*[C]{e c i} | (21) 

Ui = {eciYiCr} ) 

The expression for U 2 is not presented herein since 
U 2 is inconsequential in the development of the prop- 
erties of the continuum models. 

Step 8. The kinetic energy density of the contin- 
uum is given by 


K = ^ {uclY [m] {iid} ( 22 ) 

where a dot (•) refers to derivative with respect to 
time. 

The sign convention for the internal forces and 
generalized displacements, along with the associated 
direct stiffnesses and mass coefficients, is given in 
figures 5 and 6 for the simplified beam and plate 
continua. 

The foregoing approach is applied in the succeed- 
ing sections to the development of continuum mod- 
els for beamlike and platelike lattice trusses with pin 
joints. 


5. Application to Beamlike Lattices 


The foregoing approach is applied to the beamlike 
lattices having four longerons and an orthogonal 
tetrahedral (unsymmetric) configuration (see fig. 7). 
These trusses are primary candidates for use in the 
keel beam design of the Space Station (see fig. 8). For 
the lattice shown in figure 7, the smallest possible 
repeating cell that can be isolated from the grid 
extends over one bay of the original structure. 


5.1. Kinematic and Temperature Assumptions 


The deformed position of any cross section of the 
lattice is assumed to be specified by 12 displacement 
parameters. A Timoshenko-type beam theory is 
obtained by eliminating six of these parameters. The 
following expressions are used for the displacement 
field in the plane of the cross section (plane yz): 


{ u(x,y,z) 
v{x,y,z) 
w{x,y,z) 



where v°, and w° are the displacement com- 
ponents at y = z — 0 (chosen to be at the cen- 
ter of the repeating cell); (j) x , and (j> z are 

the rotation components about the coordinate axes; 
Sy and e° z are the extensional strains in the y- and 
^-directions; 7 ° z is the shearing strain in the plane 
of the cross section (plane yz)\ u , v , and w rep- 
resent warping and distortion of the cross section; 
and the dots refer to zero terms. The sign con- 
vention for the displacement and rotation compo- 
nents is shown in figure 5. The 12 parameters 
U°, V°, w°, 4>x, 4>y, <Pz, £y, 7 yz,u, v, and w 

are functions of x only. Note that equations (23) 
provide an exact representation of the displacement 
field in the plane of the cross section. As a conse- 
quence of the kinematic assumptions (eqs. (23)), the 
strain components have a bilinear variation in the 
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plane of the cross section (plane yz) as follows: 


£ll =e° x - yK° + zk° + yzty 0 
£22 = £y + zv 


£33 = £°z + yw 

712 = l° X y + y d x £y + z (-k? + k) + yz d x v 

713 = 7 xz+V + «) + * d x e° + yz d x w 

723 = 7 yz + yv + zw 


} (24) 


where e° is the extensional strain of the centerline; 
Ky and k° z are the curvature changes in the y- and 
2-directions; k° is the twist; 7 xy and ^ xz are the 
transverse shear strains; and d x = d/dx. The strain 
parameters £ x , k°, k°, 4zi'1xy K t > and « are 
functions of x only and can be expressed in terms of 
the displacement and rotation components as follows: 


£° = d x u° 

K y = d X <l>Z 

K z — 0 X <t>y 

7 xy = d x V° - <j> z 


7x2 = d xW° + 4>y ' 

K t = dx<t>x 
ijj° = d x u 

K = 7,d X ~iy Z + U 


(25) 


The temperature distribution is also assumed to be 
bilinear in the plane of the cross section; that is, 


T(x, y,z) = T° + y d y T° + 2 d z T° + yz d y d z T° (26) 


where T° is the temperature at y = 2 = 0 and the 
temperature gradients are 


dyT° 


dr° 

dy 


d z T° = 


<9T° 

dz 


dy d z T° = 


d 2 T° 
dy dz 


5.2. Transformation Matrices and Simplified 
Beam Models 

The kinematic assumptions (eqs. (23)) are used 

(AO 

to generate the transformation matrix of the 

fcth member. Equations (24) and (26) are used 

(AO 

to generate the transformation matrices [1^ ] and 

[rPl with y and z in these equations denoting the 
coordinates of the center of the kth member. 

A simplified (Timoshenko-type) continuum beam 
model is obtained by retaining six displacement pa- 
rameters in equations (23) and six strain parameters 
in equations (24). This is accomplished by defining 
the partitions of {s} c and {u} c as follows: 


{£cl}* = K Ky K° 
{£ c2 } t = [a 7° ; 2 tp° 
{£c3}‘ = [dx£°y d x £° z 
{UclY = [u° V° w° 
{l^c 2 } 7 yz £ 3 / 


7 °xy Ixz rf] 

£° y £° z V W] 

dxlyz d xV d x w] > 
<Px <t>y 4>z] 

£° V W] 


The constitutive matrix and the matrix of mate- 
rial densities are obtained by following the procedure 
described in section 4. The numerical values for the 
coefficients of these matrices are given in table 2 for 
the lattice structure shown in figure 7. Stiffness co- 
efficients and density parameters for the simplified 
and higher order continuum models are presented in 
reference 27 for orthogonal tetrahedral lattices with 
more than four longerons. 


6. Application to Platelike Lattices 

In this section applications of the foregoing ap- 
proach to the tetrahedral and hexahedral grids shown 
in figures 1 and 9, respectively, are outlined. The 
tetrahedral grid is chosen because it has many attrac- 
tive features for application to space structures. The 
hexahedral grid is chosen because it contains more 
redundancies and has vertical members that present 
difficulties in some continuum modeling approaches. 
In each case, a typical repeating element is isolated 
from the grid. The continuum model in this case is 
a two-dimensional plate continuum. 


6.1. Kinematic and Temperature Assumptions 

The three displacement components u, and w 
are assumed to have a linear variation in the thick- 
ness coordinate 2; that is, 


{ u(x,y,z) ) ( u° ) ( <t> x ) 

v{x, y.z) > = < v° \+zl<t>y\ (28) 

w{x,y,z)) ( J l £° J 

where u ° , v°, and w° are the displacement com- 
ponents at z = 0 (chosen to be at the center 
of the repeating cell); (j) x and (j) y are the rota- 
tion components; and e° z is the transverse nor- 
mal strain in the ^-direction. The six parameters 
u ° , v° , w°, 0 X , (fry, and e° z are functions of x and 
y only. The sign convention for the displacement 
and rotation components is shown in figure 6. Equa- 
tions (28) represent the exact displacement variation 
in the thickness direction provided that no internal 
nodes are present. 

As a consequence of the displacement assump- 
tions (eqs. (28)), the strain components have a linear 
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variation across the thickness of the plate as follows: 


' £11 1 


f £° 1 


( v-O 

K' X 

£22 


£° 

£ e l 

c z 
7 xy 


K° 

£33 

712 

► = < 

► + z < 

2 K xy 

d x e° 

713 


7 xz 


- 723 > 


7j/z J 


d y e° 2 J 


(29) 


where e ° y , and ^ xy are the extensional and shear- 
ing strains of the middle plane; k x , and 2k xv are 
the curvature changes and twist; ^ xz and 7 ^ are the 
transverse shear strains; and a dot represents zero. 
The strain measures s x ^ s y , 7 xyi ^zi^yi xy 5 'Ixzi 
and 7 y Z are independent of z and can be expressed in 
terms of the displacement and rotation components 
of the middle plane as follows: 


e° x = d x u° > 

e° y = d y v° 

7 Sy = d x v° + d v u° 

K X = dxfix 

k y — dy&y 

‘ZK'xy ~ dx4>y + dy<t>x 
1 xz = 3 *™° + <t>x 
lyz = * 


(30) 


To account for the local transverse deformation of 
the repeating cell for the platelike lattices considered 
in the present study, the transverse shear strain 
expressions must be augmented by their derivatives 
with respect to x and y as follows: 


713 = 7 xz + z d x£° z + x d x i° 2 + y d v ^° xz 1 
723 = lyz + * + * dxlyz + y d ylyz J 


The temperature distribution is also assumed to 
be linear in the ^-direction; that is, 


T(x,y,z) = T° + zd z T° (32) 

where T° is the temperature at z = 0, and 
d z T° = dT°/dz is the temperature gradient in the 
2 -direction. Both T and d z T° are functions of x and 
V Only. 


6.2. Transformation Matrices and Simplified 
Plate Models 

Equations (28), (29), (31), and (32) are used to 
generate the three transformation matrices [T^], 
[rf } ], and [T^] for the fcth member. A simplified 


shear-deformation-type plate theory is obtained by 
defining the partitions of {s c } and {u c } as follows; 

{^cl} = [^2 £ y Ixy K x K y 2 K xy ~f xz ~1yz\ 
{^c2}* = [dxlxz dylxz dxlyz dylyz £ z ] 

{^ 03 } = [dxEz > 

{UclY =[u° v° W° <t>x <t>y\ 

{tic2> t = [ez] 

(33) 

The constitutive matrix and the matrix of ma- 
terial densities are obtained by following the proce- 
dure outlined in section 4. The numerical values for 
the coefficients of these matrices are given in table 3 
for the double-layered hexahedral lattice configura- 
tion shown in figure 9. Analytical expressions for the 
stiffness coefficients and density parameters of sim- 
plified plate models are given in reference 25 for the 
double-layered tetrahedral and hexahedral grids. 

7. Comments on the Foregoing Approach 
for Developing Continuum Models 

The following comments regarding the procedure 
for developing continuum models seem to be in order: 

1. The strain parameters {£ c i}, {s C 2}, and {^ 3 } 
can be selected on the basis of a sensitivity anal- 
ysis of the response with respect to each of these 
parameters. 

2. For simple lattice configurations, the constitu- 
tive matrix [C], the matrix of density parameters [m], 
and the vector of thermal loads {Ct} can be obtained 
in symbolic form using a computerized symbolic ma- 
nipulation language such as MACSYMA (see ref. 47). 
However, for more complicated configurations, it is 
more convenient to generate [C], [m], and {C?} nu- 
merically using the FORTRAN program listed in the 
appendix on a microcomputer. 

3. Higher order continuum models can be ob- 
tained by using more strain and displacement pa- 
rameters than {e c \} and {u c 1 }. The additional pa- 
rameters may be needed for the accurate prediction 
of the detailed displacement and stress distribution 
in the lattice (for example, localized displacements 
and stresses near a loaded edge). 

4. A similar procedure for generating the char- 
acteristics of the continuum model was presented in 
references 23 and 24. In the cited references the lat- 
tice structure is considered to be the sum of several 
arrays of parallel elements, and thermoelastic stiff- 
ness coefficients of the continuum model are obtained 
by applying the tensor transformation relations to 
the unidirectional stiffness and thermal coefficients 
of each of the individual members constituting the 
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repeating cell. The three-dimensional continuum co- 
efficients generated are equal to the stiffness ma- 
trix and thermal load vector [K\ and { P T } of equa- 
tions (4) and (5), each divided by the volume of the 
repeating cell. 

8. Numerical Studies 

To test and evaluate the accuracy of the predic- 
tions of the continuum models developed, a num- 
ber of thermoelastic stress analysis, free vibration, 
and bifurcation buckling problems were solved using 
these models. Comparisons were made with exact so- 
lutions based on direct analysis of the actual lattice 
structure. Two problem sets are presented herein: 
(1) stress and free vibration analysis of anisotropic 
beamlike lattices and (2) free vibration analysis of 
platelike lattices. 

8.1. Stress and Free Vibration Analysis of 
Anisotropic Beamlike Lattices 

The first problem set considered is that of the 
orthogonal tetrahedral truss configuration shown in 
figure 7, which is a candidate for the primary truss 
support structure for the Space Station (see fig. 8). 
Because of the unsymmetry of the basic configura- 
tion, the continuum model is anisotropic in the sense 
that both the extensional and the shear effects, as 
well as the bending and twisting effects, are coupled. 
The elastic and dynamic characteristics of the contin- 
uum model are listed in table 2. For stress analysis 
problems, the structure was subjected to longitudinal 
loadings, transverse loadings, and twisting moments 
at its free end. Typical results are shown in figures 10 
and 11 and in tables 4 and 5. 

Table 4 gives the maximum displacements u ° , v ° , 
and w° at the free end obtained by the beam mod- 
els, along with the exact solutions obtained by the 
direct analysis of the actual structure for lattices 
with 10 bays. For the sake of comparison, the pre- 
dictions of the foregoing continuum model are given 
along with those of the continuum model in which 
all the coupling stiffness and density parameters are 
neglected. Table 5 and figure 10 give an indication 
of the accuracy of the lowest 10 frequencies obtained 
by the 2 continuum models. Then vibration mode 
shapes associated with the lowest 10 frequencies of 
the 10-bay truss are shown in figure 11. As can be 
seen from table 5 and figure 10, the predictions of 
the foregoing continuum models are highly accurate. 
The maximum error in the lowest 10 frequencies for 
the 10-bay truss is 4.73 percent. The error reduces 
to less than 1.5 percent when the number of bays 
increases to 30. On the other hand, neglecting the 
coupling terms results in overestimating the stiffness 


of the structure, thereby increasing the vibration fre- 
quencies and reducing the maximum displacements. 
This effect is very pronounced for the vibration fre- 
quencies associated with the higher modes. For the 
10-bay truss, when the coupling terms are neglected, 
the maximum errors in the lowest 10 frequencies are 
over 65 percent. 

8.2. Free Vibration Analysis of Platelike 
Lattice Grids 

The second problem considered is that of the 
free vibration analysis of a cantilevered hexahedral 
double-layered grid. The characteristics of the grid 
are shown in figure 9. The boundary nodes at x = 0 
are completely restrained. The continuum model is 
taken to be a square, shear-flexible plate with one 
edge completely restrained. The elastic and dynamic 
characteristics of the plate are listed in table 3. The 
continuum solutions presented herein are converged 
finite element solutions. Typical results are shown in 
figures 12 and 13 and in table 6. 

An indication of the accuracy of the lowest 10 vi- 
bration frequencies is given in table 6 and figure 12. 
The standard of comparison is taken to be the vibra- 
tion frequencies obtained by the direct finite element 
solutions of the actual grid. Figure 13 shows the vi- 
bration mode shapes associated with the minimum 
10 frequencies of the 8 L x 8 L grid. As can be seen 
from table 6, for an 8 L x 8 L grid, the maximum error 
in the lowest 10 vibration frequencies is 1.1 percent. 
The error reduces to less than 0.55 percent when a 
16 L x 16L grid is used. 

9. Application of Continuum Models to 
Stability Problems 

An important criterion in the design of light- 
weight lattice trusses is the onset of elastic instability. 
For finely divided lattice trusses the failure mode 
is usually that of general (system) instability rather 
than local snap-through of a single joint or member 
buckling. The foregoing continuum models can be 
extended to the general (global) stability analysis of 
lattice structures. This is accomplished by equating 
the potential energy due to initial stresses of the 
original lattice to that of the continuum, from which 
expressions for the geometric stiffness coefficients of 
the continuum are obtained in terms of the geometric 
and elastic properties of the original lattice structure. 
The procedure is outlined as follows: 

1. A repeating cell is isolated from the lattice 
grid. The geometric stiffness matrix of a typical 
member k of the repeating cell referred to the co- 
ordinates of the repeating element is given by [g(*)]. 
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2. The nodal displacements of member k are 
expressed in terms of the continuum displacement 
parameters and their derivatives with respect to 
the spatial coordinates by means of the following 
transformation: 


M W = [fi fc) ] {A} c (34) 

where {A} c is the vector of displacement parameters 
and their spatial derivatives and |f^J is a trans- 
formation matrix. Equations (34) are obtained by 
expanding the nodal displacements of member k in a 
Taylor series about the center of the repeating cell. 

3. The geometric stiffness matrix of the contin- 
uum is given by 


I&-S £ 

Members 



(35) 


4. Simplified continuum models are obtained by 
partitioning the vector {A} c and the matrix [G] c as 
follows: 



(36) 


I9\c = 


9 ii 9n 

_ i ?21 ^22 


(37) 


where {A c i} and {A C 2 } are the displacement param- 
eters associated with the retained and neglected geo- 
metric stiffness coefficients, respectively. The geo- 
metric stiffness coefficients of the simplified contin- 
uum are given by 


beamlike lattices. This approach is applied to static, 
large-deflection analysis and large-amplitude free vi- 
bration analysis, as well as buckling and postbuckling 
analysis. 

10.2. Beamlike Lattices With Material 
Damping 

Planar beamlike lattices with members having 
viscous damping are modeled as continuum Timo- 
shenko beams in references 52 and 53. The global 
damping characteristics are determined in terms of 
the damping coefficients and dimensions of the truss 
members. 

10.3. Stress, Free Vibration, and Buckling 
Problems of Beamlike Lattices With Rigid 
and Flexible Joints 

Micropolar continuum models are developed in 
references 48, 54, and 55 for the static, free vibra- 
tion, and buckling problems of planar and spatial 
beamlike lattices with rigid joints. The study of ref- 
erence 48 shows that ordinary shear-flexible contin- 
uum beam models are generally adequate for predict- 
ing the global response characteristics of lattices that 
do not need to have rigid joints for their kinematic 
stability. Such continuum beam models are inad- 
equate however when diagonal members have very 
small cross-sectional areas. Incorporation of joint 
flexibility in the continuum model is discussed in 
reference 56. 

10.4. Stress and Free Vibration Problems of 
Single-Layered Grids 


[?] = [9u] (38) 

Note that the foregoing continuum approach for 
stability analysis assumes that the individual mem- 
bers of the lattice remain straight and stable dur- 
ing buckling. Local member instability and snap- 
through buckling of a single joint are not predicted 
by the present theory. 


Ordinary continuum plate models are developed 
in references 20 and 57 for predicting the response of 
single-layered grids with rigid joints. The transverse 
motion of the grids is predicted accurately by these 
models. However, accurate prediction of the in-plane 
motion of single-layered grids with rigid joints re- 
quires the use of a micropolar continuum (see refs. 54, 
58, and 59). 


10. Other Reported Applications of 
Continuum Modeling 

Other reported applications of continuum model- 
ing which have not been discussed in the preceding 
sections are described subsequently. 

10.1. Geometrically Nonlinear Problems of 

Beamlike Lattices With Pin Joints 

In references 50 and 51 an incremental stiffness 
matrix of the equivalent continuum is used to ac- 
count for the large-displacement nonlinear effects in 


10.5. System Identification and Control Law 
Design of Lattice Structures 

Because of the significant reduction in the num- 
ber of material and structural parameters obtained 
by replacing the actual lattice structure by a con- 
tinuum model, a number of recent studies have em- 
ployed continuum modeling in control design and/or 
system identification analysis. These studies used ei- 
ther the governing partial differential equations of 
the continuum or the systems of ordinary differential 
equations resulting from discretizations of the con- 
tinuum model (see refs. 60 to 63). 
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11. Future Directions of Research on 
Continuum Modeling 

Continuum modeling is likely to play a significant 
role in the analysis, design, and control of large space 
structures, and accordingly, a number of extensions 
and computational strategies must be developed. 

11.1. Needed Extensions and Future 

Applications 

The needed extensions to the continuum model 
include incorporating effects of joint flexibility and 
damping. 

New applications of the continuum model include: 

1. Random periodic structures (i.e., 
structures with random variations 
in material properties, geometric 
characteristics, or periodicity) 

2. Nonlinear structural dynamics and wave 
propagation problems 

3. Multidisciplinary optimization problems 
(e.g., simultaneous optimization of 
structures and control systems) 

11.2. Computational Strategies and Modeling 
Techniques 

With regard to computational strategies and 
modeling techniques, the following two areas appear 
to have high potential for research. 

Hybrid modeling approaches. Hybrid modeling 
approaches include combination of either (1) contin- 
uum and repetitive structure approaches for handling 
lattice structures with arbitrary boundary conditions 
or (2) continuum and discrete structure approaches 
for specialized problems such as stress concentration. 

Application of operator splitting/reduction tech- 
nique for analyzing complicated continuum models . 
Operator splitting/reduction techniques are very ef- 
fective for generating the response of a complex struc- 
ture (e.g., anisotropic continuum model) using large 
perturbations from the response of a simpler struc- 
ture (e.g., corresponding orthotropic model). These 
techniques are applied to eigenvalue and nonlinear 
problems in references 64 and 65. 

Reduction techniques have been shown to be re- 
lated to the preconditioned conjugate gradient iter- 
ative method. The preconditioning matrix is taken 
to be the matrix of the simpler structure; and the 
preconditioned residuals provide sensitivity informa- 
tion of the response with respect to the complicating 
factors (e.g., anisotropy). 


12. Concluding Remarks 

The status and some recent developments in con- 
tinuum modeling of large repetitive lattice structures 
are summarized. Discussion focuses on a number of 
aspects including definition of an effective substitute 
continuum, characterization of the continuum model, 
and the different approaches for generating the con- 
stitutive matrix and the matrix of mass densities of 
the continuum. Also, a simple approach is presented 
for generating the continuum properties. The ap- 
proach can be used in conjunction with computerized 
symbolic manipulation to generate analytic expres- 
sions for the continuum properties or with a small 
FORTRAN program to determine the numerical val- 
ues of these properties. The proposed approach is 
applied to the generation of continuum properties for 
beamlike and platelike lattices. 

The status of continuum modeling for repetitive 
lattice structures can be summarized in the following: 

1. For pin-jointed beamlike and platelike lattices, 
effective and verified ordinary continuum models ex- 
ist for linear thermoelastic stress and free vibration 
problems. Bifurcation buckling loads associated with 
general instability (global buckling modes) can also 
be predicted accurately by continuum models. Ap- 
plications of continuum modeling to nonlinear static 
and dynamic problems are limited. Simplified con- 
tinuum models (associated with known engineering 
theories such as Timoshenko-type beam theory and 
Reissner-Mindlin-type plate theory) are adequate for 
many practical problems. Continuum models associ- 
ated with higher order theories may be required for 
predicting localized effects. 

2. For lattices with rigid joints, ordinary contin- 
uum models are available for the accurate prediction 
of the global response of lattices with rigid joints, 
provided that the response is not dominated by lo- 
cal member deformation (e.g., transverse motion of 
single-layered grids). Micropolar continuum models 
exist for predicting the in-plane motions of lattices 
with rigid joints. 

Needed extensions of continuum modeling in- 
clude incorporation of joint flexibility and damping. 
Also, application of continuum modeling to multidis- 
ciplinary optimization problems needs more devel- 
opment. These extensions are recommended as fu- 
ture directions for research. New computational and 
modeling strategies that have high potential for use 
in conjunction with continuum modeling include hy- 
brid combination continuum and repetitive structure 
approaches; hybrid combination of continuum and 
discrete structure approaches; and the application of 
operator splitting/reduction technique for analyzing 
complicated continuum models. 
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Appendix 

FORTRAN Program for Evaluating the Thermoelastic Stiffnesses and Density 
Parameters of the Continuum Models 

This appendix contains the listing of the interactive program used for evaluating the thermoelastic 
stiffnesses and density parameters of simplified one- and two-dimensional continuum models of beamlike 
and platelike lattice structures with pin joints. The program is written in Microsoft FORTRAN 3.30 
(1985) and works on IBM PC and XT with DOS 3.0 operating system. A tree structure depicting the 
calling sequence of the different subroutines is given in figure Al. 


CMODEL 


initial 

— menu 

mssge 

— goset — 

default 

— eqsbeam 

main 

— eqsplat 

saybye 

— fetch 


— saverun 


— run 


edsetl 

edset2 

edset3 

edset4 

edset5 


setup 
rtmgam 
kcmcpt 
partit 
cct — 



Figure Al. Tree structure for program. 


dogaml 

dogam2 

outbeam 

outplat 


li 



General Comments: 


1. There are 5 separate "include" files: 

DBLE, PARAM.CMN, DATA . CMN , CKCOLOR.CMN, GLOBCONS.CMN 

2. The following utility routines are used (the routines in 
part a were obtained from MicroCompatibles , Inc., 

301 Prelude Drive, Silver Spring, MD 20901): 

a) Grafmatic library routines: 

QBORD - set color of border 
QCLEAR - clear screen and set color 

QCMOV - move cursor to specified column and row 

b) Other routines for which the source code is not included: 

CKFIELD - read in value at specified column and row 

CKTABL - read in table of values 

CLFIELD - clear a field 

ERRORS - print error message and beep 

INITIAL - initialize key definitions and colors 

JCLEAR - clear screen and move cursor to top left corner 

JLABEL - print label at specified column and row 

PRFIELD - print value at specified column and row 

PRTABL - print table of values 

VTOASC - convert value to ascii 


« file DBLE follows» 

IMPLICIT REAL*8 (A-H, 0-Z) 

« file PARAM.CMN follows» 

C*C0MDECK PARAM 
C 

INTEGER ZNN , ZNE, ZNPT, ZNSDF , ZNDDF , ZNCN 

PARAMETER (ZNN=18, ZNE=60, ZNPT=8, ZNSDF=18, ZNDDF=12, ZNCN=2) 
C 

C ZNN = MAX. NUMBER OF NODES IN BAY 

C ZNE = MAX. NUMBER OF CONNECTING ELEMENTS IN BAY 

C ZNPT = MAX. NUMBER OF PROPERTY TYPES 

C ZNSDF = MAX. NUMBER OF STRAIN DEGREES OF FREEDOM 

C ZNDDF = MAX. NUMBER OF CONTINUUM DISPLACEMENT DEGREES OF FREEDOM 

C ZNCN = MAX. NUMBER OF NODES CONNECTED 

C 

<< file DATA. CMN follows» 

C*COMDECK CONN 
C 

COMMON /CONN/ NOP (ZNE, ZNCN) 

C 

C NOP - CONNECTIVITY TABLE, NE X NCN 
C 

C*COMDECK COORD 
C - 

COMMON /COORD/ XCOORD(ZNN) , YCOORD(ZNN) , ZCOORD(ZNN) 

C 
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C YCOORD - COORDINATES IN Y DIRECTION 
C ZCOORD - COORDINATES IN Z DIRECTION 
C 

C*COMDECK DIRCOS 
C 

COMMON /DIRCOS/ DCOSX(ZNE) , DCOSY(ZNE) , DCOSZ(ZNE) 

C 

C DCOSX - DIRECTION COSINE IN X DIRECTION (=L IN EQUATIONS) 

C DCOSY - DIRECTION COSINE IN Y DIRECTION (=M IN EQUATIONS) 

C DCOSZ - DIRECTION COSINE IN Z DIRECTION (=N IN EQUATIONS) 

C 

C*COMDECK ECUC 

C 

INTEGER EC. UC. ECI, UCI 

COMMON /ECUC/ EC(ZNSDF). UC(ZNDDF), ECI(ZNSDF), UCI(ZNDDF), 
NEC1, NEC2, NEC3, NUC1, NUC2 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


NOTE: PARAMETERS HAVE DEFAULT ORDERING AND ARE NUMBERED AS SUCH 
EC, UC, ECI. UCI TAKE CARE OF PERMUTATION OF THE PARAMETERS 


EC 


ECI 

UC 


UCI 


= (ECI, EC2, EC3) , WHERE ECI, EC2, EC3 CORRESPOND TO RETAINED. 

ASSOCIATED W/NEGLECTED, AND NEGLECTED 
PARAMETERS, RESPECTIVELY 

EC(K) = N, THE NTH PARAMETER IS IN KTH POSITION 
(INVERSE OF EC) ECI(N) = K, " " " " 

= (UCI, UC2) , WHERE UC1.UC2 CORRESPOND TO RETAINED AND 

NEGLECTED PARAMETERS, RESPECTIVELY 
UC(K) = N, THE NTH PARAMETER IS IN KTH POSITION 
(INVERSE OF UC) UCI(N) = K, " " " " 


C NEC1 , NEC2 , NEC3 , NUC1 , NUC2 - NUMBER OF PARAMETERS IN EACH GROUP 

C 

C 

C*COMDECK ECOORD 
C 

COMMON /ECOORD/ EX(ZNE) , EY(ZNE) , EZ(ZNE) 

C 

C EX - X COORDINATE AT CENTER OF ELEMENT 

C EY - Y COORDINATE AT CENTER OF ELEMENT 

C EZ - Z COORDINATE AT CENTER OF ELEMENT 

C 

C*COMDECK ELEQS 

C 

REAL*8 M 

COMMON /ELEQS/ K, R(6) , T(4), M(6,6), GAM1 (6 , ZNSDF) , GAM2 (6 , ZNDDF) 
C 

C K - WHICH NUMBER ELEMENT WORKING ON 

C [R] , [T] , [M] , [GAM1] , [GAM2] - ELEMENTAL EQUATIONS 

C 

C*COMDECK FULLEQS 
C 

PARAMETER (NSS=ZNSDF*ZNSDF , NS4=ZNSDF*4, NDD=ZNDDF* ZNDDF) 

REAL*8 KC, MC, KC11, KC12, KC21, KC22, MC11 
COMMON /FULLEQS/ KC(NSS), PT(NS4) , MC(NDD) , C(NSS), CT(NS4) , 
KCll(NSS), KC12(NSS) , 


KC21CNSS), KC22(NSS), MCll(NDD) , 

PT1(NS4), PT2(NS4) 

C 

C*COMDECK LENS 

C 

REAL*8 LEN 

COMMON /LENS/ LEN(ZNE) 

C 

C LEN - LENGTH OF EACH ELEMENT 
C 

C*COMDECK NUMBER 
C 

COMMON /NUMBER/ I TRUSS, OMEGA, NN, NE, NSDF , NDDF, NPT, NCN, NDIMT 
CHARACTER*80 TITLE 
COMMON /TITLEC/ TITLE 
C 

C TITLE - TITLE OF DATA 
C ITRUSS -=1, BEAMLIKE LATTICE TRUSS 
C =2, PLATELIKE LATTICE TRUSS 

C OMEGA - =LENGTH OF REPEATING CELL FOR BEAMLIKE LATTICE (ITRUSS=1) 

C =PLANFORM AREA OF REPEATING CELL FOR PLATELIKE LATTICE (ITRUSS=2) 

C NN - NUMBER OF NODES 

C NE - NUMBER OF ELEMENTS 

C NSDF - NUMBER OF STRAIN DEGREES OF FREEDOM 

C NDDF - NUMBER OF DISPL. DEGREES OF FREEDOM 

C NPT - NUMBER OF PROPERTY TYPES 

C NCN - NUMBER OF NODES CONNECTED IN EACH ELEMENT 

C NDIMT - DIMENSION OF ELEMENTAL MATRIX [T] , EITHER=4 OR 2 

C 

C*COMDECK PROP 
C 

COMMON /PROP/ AREA(ZNPT) , E(ZNPT) , RHO(ZNPT) , ALPHA (ZNPT) , 

IPROP(ZNE) 

C 

C AREA - CROSS-SECTIONAL AREA OF STRUCTURE 
C E - YOUNG’S MODULUS 

C RHO - MASS DENSITY 

C ALPHA - COEFFICIENT OF THERMAL EXPANSION 
C IPROP - FOR EACH ELEMENT, WHICH PROPERTY TYPE 
C 

« File GLOBCONS . CMN follows » 

common /globcons/ SOLID, DASHED, DOTTED, LEFT, RIGHT, MIDDLE, 

* ARROW, CROSS, BLOCK, POINT, HRGLAS, 

* YES, NO, ON, OFF, TRUE, FALSE, 

* BLACK. BLUE. GREEN. CYAN, RED, MAGENT, 

* BROWN. LGRAY, DGRAY, LBLUE , LGREEN , LCYAN . 

* LRED, LMAGEN , YELLOW, WHITE, ULINE 
integer SOLID, DASHED. DOTTED, RIGHT, ARROW, CROSS, BLOCK, 

* POINT, HRGLAS, YES, ON, OFF, TRUE, FALSE, 

* BLACK, BLUE, GREEN, CYAN. RED, BROWN, DGRAY, YELLOW, 

* WHITE, ULINE 


« file CKCOLOR . CMN follows » 
c 
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c BG=BACKGROUND , FG=FOREGROUND , T=TEXT , F=FIELD, W=WINDOW, 
c IC0L0R=1 FOR COLOR MONITOR, 0 FOR MONOCHROME 
c 


* 


common /ckcolor/ IBGT, IFGT, IBGF , IFGF, IBGTW , IFGTW , 
IBGFW , IFGFW, ICOLOR 


oooooooooooooonoo 


$ STORAGE : 2 

C*******************************^************************************ 
C CONTINUUM MODEL FOR REPETITIVE LATTICE STRUCTURES: 

C 

C BEAMLIKE LATTICE TRUSSES 

C PLATELIKE LATTICE TRUSSES 

C 

*** USES DOUBLE PRECISION *** 

GENERAL EQUATIONS: 


FOR EACH ELEMENT - 

[R] = [L*L, M*M, N*N, L*M, L*N, M*N] 

[1 0 0] 

[M] = 1/6 * RHO*AREA*LEN * [2*13 13] 13= [0 1 0] 

[ 13 2*13] , [001] 

[T] = [1, YC, ZC, YC*ZC] , BEAMLIKE 

[T] = [1, ZC] , PLATELIKE 

[GAM1] , DEPENDS ON BEAM-OR-PLATE, CHOICE OF 
[GAM2] , EC1 , EC2 , EC3 , UC1 , UC2 
C 
C 

C [KC] = SUM OF ELEMENTS ( [R] * [GAM1] )T * ([R]*[GAM1] 

C [PT] = SUM OF ELEMENTS ( [R] * [GAM1] )T * [T] 

C [MC] = SUM OF ELEMENTS 1/OMEGA * [GAM2] T* [M] * [GAM2] 

C 

C [KC] = [ KC11 KC12 KC13 ] [MC] = [MC11 MC12] 

C [ KC22 KC23 ] [MC21 MC22] 

C [ SYM. KC33 ] 

C 

C [C] = 1/OMEGA * [KC11 - (KC12*KC22 -1 ) * KC21] 

C 

C [CT] = 1/OMEGA * [PT1] - [KC12*KC22 _1 ] * [PT2] 

C 

C 

C FILES : 

C 

C UNIT 5 - SCREEN INPUT 

C UNIT 6 - PROBLEM OUTPUT 

C UNIT 7 - DEBUG FILE (NO LONGER USED) 

C UNIT 8 - FILE TO READ DATA (OR TO WRITE DATA AND RESULTS) 

Q****** ********** ********* ************************************** ******** 

SINCLUDE: ’DBLE’ 

C unit 5 - screen 

C unit 6 - screen 

OPEN (5,FILE=’C0N: ’) 

OPEN (6,FILE=’C0N: ’) 

c initialize key definitions and colors . 
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CALL INITIAL 


c . write cover message . 

CALL MSSGE 

c . initial variables and set defaults . 
CALL DEFAULT 

c . start main job . 

CALL MAIN 

c . write ending message . 

CALL SAYBYE 
END 


SUBROUTINE MSSGE 

c * * * * * *** * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * ** 
C Writes cover message and any instructions 

C * * ** **** * *** ** * * * *** ****** *********** *** *** ****** *** ********** ****** * * * 

$INCLUDE: ’DBLE’ 

$ INCLUDE : ’ GLOBCONS . CMN ’ 

$ INCLUDE: ’ CKCOLOR . CMN ’ 

CALL QCLEAR (BLUE , YELLOW) 

CALL QBORD (GREEN) 

CALL QCMOV (0,24) 


CALL JLABEL (12. 20, 
. ctures $’) 

CALL JLABEL (12, 19, 

$’) 

CALL JLABEL (5. 17, 
.owing matrices : $’) 
CALL JLABEL (15. 15. 
. [C] $') 

CALL JLABEL (15, 14, 
. . [CT] $’) 

CALL JLABEL (15, 13, 
.ontinuum, [m]$’) 
CALL JLABEL ( 5, 11, 
CALL JLABEL (15, 9, 

CALL JLABEL (15, 8, 


’Continuum Models for Repetitive Lattice Stru 


’Welcome ! This program will compute the foil 

stiffness matrix of equivalent continuum, 

’thermal coefficients of equivalent continuum 

’matrix of density parameters of equivalent c 

’for the following structures : $’) 

’beamlike lattice trusses $’) 

’platelike lattice trusses $’) 


CALL QCMOV (0,2) 

PAUSE 

RETURN 

END 


SUBROUTINE DEFAULT 

C Sets defaults (or do initializing) for data 

Q* ********************************************************************** 
$ INCLUDE: ’DBLE’ 

$ INCLUDE: ’ PARAM . CMN ’ 



o o o 


SINCLUDE: ’DATA.CMN’ 


C — data set 1 — 

ITRUSS = -1 
OMEGA = -1 
NN = -1 
NE = -1 
NPT = -1 
NCN = 2 
NSDF = 0 
NDDF = 0 

C — data set 2 — 

DO 10 I = 1, ZNE 
DO 10 J = 1, ZNCN 
NOP(I.J) = 0 
10 CONTINUE 

DO 20 I = 1, ZNE 
IPROP(I) = 1 
20 CONTINUE 

DO 30 I = 1, ZNPT 
AREA(I) = 0. 

E(I) = 0. 

RHO(I) = 0. 

ALPHA ( I) = 0. 

30 CONTINUE 

C — data set 3 — 

NEC1 = 0 
NEC2 = 0 
NEC3 = 0 
NUC1 = 0 
NUC2 = 0 

DO 40 I = 1. ZNSDF 
EC(I) = I 
40 CONTINUE 

DO 50 I = 1, ZNDDF 
UC(I) = I 
50 CONTINUE 

C — data set 4 — 

DO 60 I = 1, ZNN 
XCOORD(I) = 0. 

YCOORD(I) = 0. 

ZCOORD(I) = 0. 

60 CONTINUE 
RETURN 
END 

SUBROUTINE MAIN 

0 **************************************************************** 


Main subroutine to control outer menu and response 
q********** ********************************* ********************* 
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$ INCLUDE: ’ DBLE ’ 

$ INCLUDE: ’GLOBCONS.CMN’ 

$ INCLUDE: ’ CKCOLOR . CMN ’ 

CHARACTER CH0ICE*5, YN*1, FIRST*1 
INTEGER INUM(l) 

COMMON /LINES/ LINES (10) 

DATA CHOICE / ’ ’/ 

1 CALL MENU 

2 CALL CLFIELD (9.24.71.IBGF, IFGF) 

CALL JLABEL (9,24, ’What do you want to do ? $’) 

CALL CKFIELD (34, 24, 5, CHOICE, ’CH’) 

C convert 1st character to ascii to see if number 
FIRST = CH0ICE(1 : 1) 

CALL VTOASC (INUM(l), 1. ’CH’, FIRST) 

INUM(l) = INUM(l) - 48 

IF (INUM(l) . GE . 1 .AND. INUM(l) . LE. 5) THEN 
CALL GOSET(INUM(D) 

GOTO 1 

ELSE IF (CHOICE ,EQ. ’B’ .OR. CHOICE .EQ. ’b’) THEN 
CALL EQSBEAM 
GOTO 1 

ELSE IF (CHOICE .EQ. ’P’ .OR. CHOICE .EQ. ’p’) THEN 
CALL EQSPLAT 
GOTO 1 

ELSE IF (CHOICE .EQ. ’F’ .OR. CHOICE .EQ. ’f’) THEN 
CALL FETCH 
GOTO 1 

ELSE IF (CHOICE .EQ. ’R’ .OR. CHOICE .EQ. ’r’) THEN 
CALL RUN 
GOTO 1 

ELSE IF (CHOICE .EQ. ’S’ .OR. CHOICE .EQ. ’s’) THEN 
CALL SAVERUN 
GOTO 1 

ELSE IF (CHOICE .Eq. ’QUIT ’ .OR. CHOICE .EQ. ’quit ’) THEN 
LINES(l) = 0 

CALL CLFIELD (9, 24 , 71 , IBGF , IFGF) 

CALL JLABEL (9 ,24, ’Are you sure you want to quit (Y/N) ? $’) 
CALL CKFIELD (50, 24,1, YN.’CH’) 

IF (YN .NE. ’Y’ .AND. YN .NE. ’y’) GOTO 2 
ELSE 

LINES (1) = 0 

CALL ERRORS (’That was not a choice ! Hit enter & try again. $’) 
CALL WAIT 
GOTO 2 
END IF 
RETURN 
END 
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SUBROUTINE MENU 

C * * * * * * * * * * * * * ** * * ** * * * * * * * * ** * ** * * * * * * * * % ** * * * * * * % * * * ** * * * * * * *** 

c 

C subroutine to print the menu screen 
C 

SINCLUDE: ’ DBLE ’ 

$ INCLUDE: ’GLOBCONS.CMN’ 

$ INCLUDE: ’ CKCOLOR . CMN ’ 


CALL JCLEAR 

CALL JLABEL (0.24. ’MENU - !’) 

CALL JLABEL (0,20, ’ENTER : $’) 

CALL JLABEL (9. 19. ’EITHER number of data set to EDIT data$’) 

DO 100 IROW = 17,13,-1 

CALL CLFIELD ( 13 . IROW . 53 . BLACK , YELLOW) 

100 CONTINUE 

CALL JLABEL (15, 17, ’SET # 1 — DESCRIPTION OF TRUSS$’) 

CALL JLABEL (15, 16, ’SET # 2 — ELEMENT DATA k PROPERTY TYPES$’) 

CALL JLABEL (15, 15, ’SET # 3 — COORDINATES!’) 

CALL JLABEL (15, 14, ’SET # 4 — CONTINUUM STRAIN PARAMETERS!’) 

CALL JLABEL (15, 13, ’SET # 5 — CONTINUUM DISPL. PARAMETERS!’) 


CALL JLABEL (9 , 10 , ’ OR $ ’ ) 
CALL JLABEL (9. 8,’"b" 

.ce truss !’) 

CALL JLABEL (9. 7,”'p'' 

.ice truss !’) 

CALL JLABEL (9, 6,’"f" 
CALL JLABEL (9. 5,’"r" 
CALL JLABEL (9. 4,”'s" 
CALL JLABEL (9. 3. ’"quit" 
RETURN 
END 


- to display equations for BEAMlike latti 

- to display equations for PLATElike latt 

- to FETCH data from a file !’) 

- to RUN data !’) 

- to SAVE data and results to file !’) 

- to QUIT !’) 


SUBROUTINE WAIT 

C** ************************************************************* ******** 
C waits until <enter> is hit 

C*********************************************************************** 

CHARACTER CH*1 

CALL CKFIELD (78, 24, 1, CH, ’CH’) 

RETURN 

END 


SUBROUTINE EQSBEAM 

C displays equations assumed for beamlike lattice truss 

C******** ****** ****** *********** ********** ****************** ************ 

! INCLUDE: ’DBLE’ 

SINCLUDE: ’GLOBCONS.CMN’ 

! INCLUDE: ’ CKCOLOR . CMN ’ 
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CALL JCLEAR 

CALL JLABEL (15,22, ’***** EQUATIONS for BEAMLIKE LATTICE TRUSS ** 

, *** $ ’ ) 

CALL JLABEL (0, 20, ’Displacement Assumptions $’) 

CALL JLABEL (0, 19, ’ $’) 

CALL JLABEL (0, 18, ’u = uo - y*phiz + z*phiy + y*z*ubar $’) 

CALL JLABEL (0, 17, ’v = to + y*epsyo + z*(- phix+gammayzo/2) + y*z 
. *vbar $’) 

CALL JLABEL (0, 16, ’w = wo + z*epszo + y* (phix+gammayzo/2) + y*z 
. *wbar $’) 

CALL JLABEL (0, 14, ’Strain Assumptions $’) 

CALL JLABEL (0, 13, ’ - $’) 

CALL JLABEL (0, 12, ’epsli = epsxo - y*kappayo + z*kappazo + y*z*p 

si$’) 

CALL JLABEL (0, 11, ’eps22 = epsyo + z*vbar $’) 

CALL JLABEL (0, 10, ’eps33 = epszo + y*wbar $’) 

CALL JLABEL (0, 9, ’gammal2 = gammaxyo + y*d(epsyo)/dx + z*(-kapp 

.ato+kappabr) + y*z*d(vbar)/dx $’) 

CALL JLABEL (0, 8, ’gammal3 = gammaxzo + z*d (epszo) /dx + y*( kapp 

.ato+kappabr) + y*z*d(wbar)/dx $’) 

CALL JLABEL (0, 7, ’gamma23 = gammayzo + y*vbar + z*wbar $’) 

CALL JLABEL (0, 5, ’Temperature Assumptions $’) 

CALL JLABEL (0, 4, ’ $’) 

CALL JLABEL (0, 3, ’T = To + y*d(To)/dy + z*d(To)/dz + yz*d2(To)/ 

.dydz $’) 

CALL CLFIELD (0 , 24 , 80 , IBGF , IFGF) 

CALL JLABEL (0,24, ’Hit enter to continue ... $’) 

CALL WAIT 
CALL JCLEAR 
RETURN 
END 


SUBROUTINE EQSPLAT 

c*********************************************************************** 
C displays equations assumed for platelike lattice truss 
c* ********************************************************************** 
$ INCLUDE: ’DBLE’ 

$ INCLUDE: ’ GLOBCONS . CMN ’ 

$ INCLUDE: ’ CKCOLOR . CMN ’ 

CALL JCLEAR 

CALL JLABEL (15,22, ****** EQUATIONS for PLATELIKE LATTICE TRUSS * 

. **** $ ’ ) 

CALL JLABEL (0, 20, ’Displacement Assumptions $’) 

CALL JLABEL (0, 19, ’ — $’) 

CALL JLABEL (0, 18, ’u = uo + z*phix $’) 

CALL JLABEL (0, 17, ’v - vo + z*phiy $’) 

CALL JLABEL (0, 16, ’w » wo + z*epsz $’) 

CALL JLABEL (0, 14, 'Strain Assumptions $’) 
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CALL JLABEL (0, 13, ’ $’) 

CALL JLABEL (0, 12, ’epsll = epsxo + z*kappaxo $’) 

CALL JLABEL (0, 11, ’eps22 = epsyo + z*kappayo $’) 

CALL JLABEL (0, 10, ’eps33 = epszo $’) 

CALL JLABEL (0, 9, ’gammal2 = gammaxyo + z*2kappaxyo $’) 

CALL JLABEL (0, 8, ’gammal3 = gammaxzo + z*d(epszo)/dx + x*d(gamm 

.axzo)/dx + y*d (gammaxzo )/dy $’) 

CALL JLABEL (0, 7, ’gamma23 = gammayzo + z*d(epszo)/dy + x*d(gamm 

.ayzo)/dx + y*d (gammayzo) /dy $’) 

CALL JLABEL (0, 5, ’Temperature Assumptions $’) 

CALL JLABEL (0, 4, ’ $’) 

CALL JLABEL (0, 3, ’T = To + z*d(To)/dz $’) 

CALL CLFIELD (0 , 24 ,80 , IBGF, IFGF) 

CALL JLABEL (0,24,’ Hit enter to continue ... $’) 

CALL WAIT 
CALL 'JCLEAR 
RETURN 
END 


SUBROUTINE FETCH 

C Fetches a data set on a file 

C (data in same format as written to in ’SAVERUN’) 

$ INCLUDE: ’ DBLE ’ 

$INCLUDE: ’ GLOBCONS . CMN ’ 

$ INCLUDE : ’ CKCOLOR . CMN ’ 

$ INCLUDE: ’ PARAM . CMN ’ 

$ INCLUDE: ’DATA. CMN’ 

CHARACTER FILENAM*12 , STRNG1*38, STRNG2*8, STRNG3*17 

LOGICAL OLDFILE 

COMMON /LINES/ LINES (10) 

c — read in name of file with data — 

CALL JCLEAR 

CALL JLABEL (0,16, ’Name of the file with data : $’) 

5 CALL PRFIELD (40,16,12, ’ \ ’CH’ , IBGF. IFGF) 

CALL CKFIELD (40, 16, 12.FILENAM, ’CH’) 

INQUIRE (FILE=FILENAM, EXIST=OLDFILE, IOSTAT=IERR) 

IF (.NOT. OLDFILE) THEN 
LINES(l) = 0 

CALL ERRORS (’That file does not exist! Hit enter. $’) 
CALL WAIT 
RETURN 
END IF 

OPEN (8 ,FILE=FILENAM , FORM= ’ FORMATTED ’ , STATUS* ’ OLD ’ ) 

c — read data to file — 

READ (8. 6) STRNG1 , ITRUSS 
READ (8, 6) STRNG1 , NN 
READ (8, 6) STRNG1 , NE 
READ (8, 6) STRNG1 , NPT 
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READ (8, 7) STRNG1 , OMEGA 

READ (8, 10) (I, (KOP(I.J). J=1,NCN). IPROP(I) . 1=1, NE) 
READ(8,20) (I, AREA(I) , E(I) , RHO(I), ALPHA ( I) , 1=1, NPT) 
READ (8 , 30) STRNG2 , NEC1 
READ (8, 40) STRNG3 , (ECU), I=1,NEC1) 

READ (8. 30) STRNG2 , NEC2 

READ (8 ,40) STRNG3 , (EC(I) . I=NEC1+1 , NEC1+NEC2) 

READ (8 ,30) STRNG2 , NEC3 

READ (8, 40) STRNG3 , (ECU), I=NEC1+NEC2+1 .NSDF) 

READ (8, 30) STRNG2 , NUC1 

READ (8, 40) STRNG3 , (UC(I) , I*1.NUC1) 

READ (8, 30) STRNG2 , NUC2 

READ (8, 40) STRNG3 , (UC(I) . I=NUC1+1 , NUC1+NUC2) 

READ (8, 50) (J. XCOORD(J) , YCOORD(J) , ZCOORD(J), J=1,NN) 
CLOSE (8) 

6 FORMAT (A. 15) 

7 FORMAT (A, E17.10) 

10 FORMAT (IX/ (15, 3X. 215, 10X, 15)) 

20 FORMAT (1X/1X/ (I6,2X,4E17. 10)) 

30 FORMAT (A, 15) 

40 FORMAT (A, (2014)) 

50 FORMAT (IX/ (15, 5X, 3E17.10)) 

C — set other variables — 

IF (ITRUSS.EQ.l) THEN 
NSDF = 18 
NDDF =12 
ELSE 

NSDF = 15 
NDDF = 6 

END IF 

CALL CLFIELD (0 , 24 , 80 , IBGF , IFGF) 

CALL JLABEL (9,24, ’Data read in. Hit enter. $’) 

CALL WAIT 

RETURN 

END 


SUBROUTINE SAVERUN 

C saves data and results of run to a file 

C*********************************************************************** 

$INCLUDE: ’ DBLE ’ 

$ INCLUDE : ’ GLOBCONS . CMN ’ 

$INCLUDE : ’ CKCOLOR . CMN ’ 

SINCLUDE: ’ PARAM . CMN ’ 

$ INCLUDE: 'DATA. CMN’ 

CHARACTER FILENAM*12, ANS*3 

LOGICAL OLDFILE 

COMMON /LINES/ LINES (10) 

c — read in name to call file — 

CALL JCLEAR 

CALL JLABEL (0,16, ’Name of the file to save results to : $’) 
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5 CALL PRFIELD (40,16,12, * ’ , ’ CH \IBGF , IFGF) 

CALL CKFIELD (40 ,16,12, FILENAM , ’ CH ’ ) 

INQUIRE (FILE=FILENAM, EXIST=0LDFILE, I0STAT=IERR) 

IF (OLDFILE) THEN 
LINES(l) = 0 

CALL ERRORS (’That file already exists! Hit enter and try again 

.$’) 

CALL WAIT 
GOTO 5 
ENDIF 

OPEN (8 . FILE=FILENAM , FORM= ’ FORMATTED ’ , STATUS= ’ NEW ’ ) 

c — write data to file — 

WRITE (8, 6) ’ITRUSS (Beamlike-1 , Platelike=2) = ’, ITRUSS 

WRITE (8, 6) ’Number of nodes in repeating cell = ’ , NN 

WRITE (8, 6) ’Number of connecting elements = ’ , NE 

WRITE (8, 6) ’Number of property types = ’, NPT 

WRITE (8 ,7) ’Characteristic dimension of cell = ’, OMEGA 

WRITE (8, 10) (I, (NOP(I.J). J=1 , NCN) , IPROP(I) . 1=1, NE) 

WRITE (8, 20) (I, AREA(I) , E(I) , RHO(I) , ALPHA (I) , 1=1, NPT) 

WRITE (8, 30) NEC1 , (EC(I) , I=1,NEC1) 

WRITE(8 ,31) NEC2, (EC(I), I=NEC1+1.NEC1+NEC2) 

WRITE (8, 32) NEC3, (EC(I), I=NEC1+NEC2+1 , NSDF) 

WRITE (8, 40) NUC1 , (UC(I) , I=1,NUC1) 

WRITE (8, 41) NUC2, (UC(I) , I=NUC1+1 , NUC1+NUC2) 

WRITE (8, 50) (J, XCOORD(J) , YCOORD(J), ZCOORD(J) , J=1,NN) 

WRITE (8,*) ’ ’ 

6 FORMAT (A, 15) 

7 FORMAT (A, E17.10) 

10 FORMAT ( ’ OELEMENT NOP PROPERTY TYPE’/ 

(15, 3X, 215, 10X, 15)) 

20 FORMAT (’OPROPERTY TYPE - ’/ 

’ TYPE NO. AREA E RHO 

ALPHA’/ 

(I6,2X,4E17. 10)) 

30 FORMAT (’0NEC1 = ’ , 15/ ’ VECTOR (EC1) IS ’, (2014)) 

31 FORMAT ( ’ 0NEC2 = ’ , 15/ ’ VECTOR (EC2) IS ’. (2014)) 

32 FORMAT ( ’ 0NEC3 = ’ . 15/ ’ VECTOR (EC3) IS ’. (2014)) 

40 FORMAT (’0NUC1 = ’ , 15/ ’ VECTOR (UC1) IS ’, (2014)) 

41 FORMAT ( ’ 0NUC2 = ' , 15/ ’ VECTOR (UC2) IS ’, (2014)) 

50 FORMAT (’ONODAL COORDINATES (X,Y,Z)’/ (15. 5X, 3E17.10)) 

c — write results to file — 

CALL OUTMATX (C, NEC1, NEC1, NEC1, 10. 8, ’ [C] : ’) 

CALL OUTMATX (CT, NEC1, NEC1, NDIMT , 10, 8, ’ [CT] : ’) 

CALL OUTMATX (MC11 ,NUC1 ,NUC1 , NUC1, 10, 8, ’ [m] : ’) 

CALL CLFIELD (0 , 24 , 80 , IBGF , IFGF) 

CALL J LABEL (9,24, ’Data and results saved. Hit enter. $’) 

CALL WAIT 
CLOSE (8) 

RETURN 

END 


SUBROUTINE SAYBYE 
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c************************************************************ *********** 
C writes ending message and clears screen 

Q*********************************************************************** 

$ INCLUDE: ’DBLE’ 

SINCLUDE: ’ GLOBCONS . CMN * 

$ INCLUDE: ’ CKCOLOR . CMN ’ 

CALL QCLEAR (BLUE. YELLOW) 

CALL QBORD (GREEN) 

CALL QCMOV (0,24) 

CALL JLABEL (19. 10,’ Hope this program has been helpful ! $’) 

CALL JLABEL (28, 8, ’Have a good day ! $’) 

CALL QCMOV (0,2) 

PAUSE 

CALL QCLEAR (BLACK. WHITE) 

CALL QBORD (BLACK) 

CALL QCMOV (0,24) 

RETURN 

END 


$STORAGE : 2 

SUBROUTINE GOSET(ISET) 

C******************************************* **************************** 

c 

C Subroutine to control data set number ISET 
C 

$ INCLUDE: ’ DBLE ’ 

CHARACTER* 9 SET (5) 

SET(l) = ’SET 1 : $’ 

SET(2) = ’SET 2 : $ ’ 

SET (3) = ’SET 3 : $’ 

SET (4) = ’SET 4 : $ ’ 

SET (5) = ’SET 5 : $ ’ 

IF (ISET .LT. 1 .OR. ISET .GT. 5) THEN 
WRITEC5, *) ’ERROR ! ISET = ’.ISET 
STOP 
ENDIF 

CALL JCLEAR 

CALL JLABEL (0,24, SET (ISET)) 

GOTO (1,2, 3, 4, 5) ISET 

1 CALL EDSET1 
GOTO 90 

2 CALL EDSET2 
GOTO 90 

3 CALL EDSET3 
GOTO 90 

4 CALL EDSET4 
GOTO 90 

5 CALL EDSET5 
90 RETURN 

END 

SUBROUTINE EDSET1 

C edits data set 1 - description of truss 

Q* **************************************************************** ****** 

SINCLUDE: ’DBLE’ 

$INCLUDE: ’ GLOBCONS . CMN ’ 

$ INCLUDE: ’ CKCOLOR . CMN ’ 

$ INCLUDE: ’ PARAM . CMN ’ 

$ INCLUDE: ’DATA. CMN’ 

COMMON /LINES/ LINES (10) 

CALL JLABEL (0,18, ’1 - Beamlike truss=l, Platelike truss=2 $’) 

CALL JLABEL (0,16, ’2 - Number of nodes in typical repeating cell, 
.NN (max=18) $’) 

CALL JLABEL (0,14, ’3 - Number of connecting elements in repeating 
.cell, NE (max=60) $’) 
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CALL JLABEL (0,12, ’4 - Number of property types, NPT 
(max=8) $’) 

CALL JLABEL (0,10,’ 5 - Characteristic dimension of typical repeati 
.ng cell $’) 

CALL JLABEL (0, 9,’ (length for beamlike, planform area for plate 
.like lattice truss) $’) 

ICOL = 68 

IF (ITRUSS.EQ.-l) CALL PRFIELD (ICOL, 18, 3,’ ’, ’ CH ’ , IBGF . IFGF) 

IF (NN.EQ.-l) CALL PRFIELD (ICOL, 16. 3,’ ’CH’ , IBGF , IFGF) 

IF (NE.EQ.-l) CALL PRFIELD (ICOL, 14, 3,’ ’. ’CH’ , IBGF, IFGF) 

IF (NPT. EQ . -1) CALL PRFIELD (ICOL, 12, 3,’ ’, ’CH’ , IBGF, IFGF) 

IF (OMEGA. EQ.-l) CALL PRFIELD (ICOL. 10,10,’ \ ’CH’. 

IBGF, IFGF) 

IF (ITRUSS. NE.-l) CALL PRFIELD (ICOL, 18, 3 , ITRUSS, ’12’ . IBGF, IFGF) 

IF (NN.NE.-l) CALL PRFIELD (ICOL, 16, 3, NN ,’ 12 ’, IBGF , IFGF) 

IF (NE.NE.-l) CALL PRFIELD (ICOL, 14, 3, NE ,’ 12 ’, IBGF , IFGF) 

IF (NPT.NE.-l) CALL PRFIELD (ICOL, 12, 3, NPT, ’ 12’ , IBGF, IFGF) 

IF (OMEGA. NE. -1) CALL PRFIELD (ICOL, 10, 10, OMEGA, ’R8’ .IBGF, IFGF) 

10 CALL CKFIELD (ICOL, 18, 3, ITRUSS, ’12’) 

IF (ITRUSS. NE.l .AND. ITRUSS. NE. 2) THEN 
LINES (1) = 0 

CALL ERRORS (’Please input 1 or 2. Hit enter to redo. $’) 

CALL WAIT 
GOTO 10 
END IF 

IF ( ITRUSS. EQ.l) THEN 
NSDF = 18 
NDDF = 12 
ELSE 

NSDF = 15 
NDDF = 6 

END IF 

20 CALL CKFIELD (ICOL, 16, 3, NN, ’12’) 

IF (NN.LT.l .OR. NN.GT.ZNN) THEN 
LINES (1) = ZNN 
LINES (2) = 0 

CALL ERRORS (’Please redo: hit enter and input an integer betwe 
.en 1 and $’) 

CALL WAIT 

GOTO 20 fc 

ENDIF 1 

30 CALL CKFIELD (ICOL, 14, 3, NE, ’12’) 

IF (NE.LT.l .OR. NE.GT.ZNE) THEN 
LINES (1) = ZNE 
LINES (2) = 0 

CALL ERRORS (’Please redo: hit enter and input an integer betwe 
.en 1 and $’) 

CALL WAIT 
GOTO 30 
ENDIF 
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40 CALL CKFIELD (ICOL, 12, 3, NPT, ’12’) 

IF (NPT.LT.l .OR. NPT.GT. ZNPT) THEN 
LINES (1) = ZNPT 
LINES (2) = 0 

CALL ERRORS (’Please redo: hit enter and input an integer betwe 
.en 1 and $’) 

CALL WAIT 
GOTO 40 
ENDIF 

50 CALL CKFIELD (ICOL, 10, 10, OMEGA, ’R8’) 

RETURN 

END 


SUBROUTINE EDSET2 

C edits data set 2 - connectivities and prop. types 

SINCLUDE: ’ DBLE ’ 

$ INCLUDE: ’ GL0BC0NS . CMN ’ 

$INCLUDE: ’ CKC0L0R . CMN ’ 

$INCLUDE: ’ PARAM . CMN ’ 

$ INCLUDE: ’DATA. CMN’ 

PARAMETER (NC0LS1=3, NC0LS2=4) 

PARAMETER (NVAL1=ZNE*ZNCN+ZNE, NVAL2=ZNPT*NC0LS2) 

COMMON /LINES/ LINES (10) 

INTEGER IlCOLSP(NCOLSl) , FLP1 (NC0LS1) , VALl(NVALl) 

, I2C0LSP(NC0LS2) , FLP2(NC0LS2) 

REAL*8 VAL2(NVAL2) 

CHARACTER TYPEP1(NC0LS1)*2, TYPEP2(NC0LS2)*2 
LOGICAL HOMO, ERR1, ERR2 
DATA I1C0LSP / 18, 24, 40/ 

DATA FLP1 /3,3,3/ 

c — check that number of elements is defined — 

IF (NE.EQ.-l) THEN 
LINES (1) = 0 

CALL ERRORS (’Please define NE in data set #1 first. Hit enter 

.to con tinue •') 

CALL WAIT 

RETURN 

ENDIF 

c — print note — 

CALL JLABEL (50,11, ’ *** NOTE *** $’) 

CALL JLABEL (50,10, ’For elements shared by n $’) 

CALL JLABEL (50, 9, ’repeating cells, the cross-$’) 

CALL JLABEL (50, 8, ’sectional areas are to be $’) 

CALL JLABEL (50, 7, ’divided by n. $’) 
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— property table — 

CALL JLABEL (14, 22, ’cross sectional Young’s modulus mass dens 
.ity coefficient of thermal$’) 



Area E 


CALL JLABEL (0,21, ’Prop. type no. 

Rho expansion Alpha$’) 

DO 30 I * 1, NPT 

VAL2(I) = AREA(I) 

VAL2(I+NPT) = E(I) 

VAL2(I+2*NPT) = RHO(I) 

VAL2(I+3*NPT) = ALPHA (I) 

30 CONTINUE 

HOMO = .TRUE. 

TYPEP2(1) = ’R8’ 

I2C0LSP(i) = 18 
FLP2(1) = 14 

NROWS = NPT 
NLINES = 4 

CALL PRTABL (0. 17, NLINES, NROWS, NC0LS2 , I2C0LSP , FLP2, VAL2, 
TYPEP2 , HOMO) 

CALL CKTABL (0, 17, NLINES, NROWS, NC0LS2 , I2C0LSP , FLP2. VAL2. 

TYPEP2 , HOMO) 

DO 40 I = 1. NPT 
AREA(I) = VAL2(I) 

E(I) = VAL2(I+NPT) 

RHO (I) = VAL2(I+2*NPT) 

ALPHA ( I) = VAL2(I+3*NPT) 

40 CONTINUE 

— element table — 

CALL JLABEL (0, 14. ’Element no. Connectivities Property type 
.$’) 

DO 10 I = 1. NE 

VALl(I) = NOP(I.l) 

VALlG+NE) = N0P(I,2) 

VAL1(I+2*NE) = IPROP(I) 

10 CONTINUE 

HOMO = .FALSE. 

DO 20 I = 1. NC0LS1 
TYPEP1 (I) = ’12’ 

20 CONTINUE 
NROWS = NE 
NLINES = 10 

22 CALL PRTABL (0. 4. NLINES, NROWS. NC0LS1. I1C0LSP, FLP1, VAL1. 
TYPEP1 , HOMO) 

CALL CKTABL (0, 4, NLINES, NROWS. NC0LS1, I1C0LSP, FLP1. VAL1, 
TYPEP1 , HOMO) 

ERR1 = .FALSE. 

ERR2 = .FALSE. 

DO 25 I = 1, NE 

NOP(I.l) = VALl(I) 

NOP (1, 2) = VALKI+NE) 

IPROP(I) = VAL1 (I+2*NE) 

IF (NOP(I.l) .LT.l .OR. N0P(I,1) .GT.NN .OR. 

N0P(I,2) .LT.l .OR. NOP (I, 2) .GT.NN) ERR1=.TRUE. 

IF (IPROP(I) .LT.l .OR. IPROP(I) .GT.NPT) ERR2= . TRUE . 

25 CONTINUE 

IF (ERR1) THEN 
LINES (1) = NN 



LINES (2) = 0 

CALL ERRORS (’Please redo: connectivity numbers should be integ 
.ers between 1 and $’) 

CALL WAIT 
GOTO 22 
ENDIF 

IF (ERR2) THEN 
LINES(l) = NPT 
LINESC2) = 0 

CALL ERRORS (’Please redo: property numbers should be integers 
.between 1 and $’) 

CALL WAIT 
GOTO 22 
ENDIF 
RETURN 
END 


SUBROUTINE EDSET3 

C edits data set 3 - coordinates 

$ INCLUDE: ’ DBLE ’ 

$ INCLUDE: ’ GLOBCONS . CMN ’ 

$INCLUDE: ’ CKCOLOR . CMN ’ 

$ INCLUDE: ’ PARAM . CMN ’ 

$ INCLUDE: ’DATA. CMN’ 

PARAMETER (NC0LS=3, NVAL=NCOLS*ZNN) 

INTEGER ICOLSP(NCOLS) , FLP(NCOLS) 

CHARACTER TYPEP(NCOLS) *2 
REAL*8 VAL(NVAL) 

LOGICAL HOMO 

COMMON /LINES/ LINES (10) 

c — check that NN is defined — 

IF (NN.EQ.-l) THEN 
LINES (1) = 0 

CALL ERRORS (’Please define NN in data set #1 first. Hit enter 
.to continue $’) 

CALL WAIT 
RETURN 
ENDIF 


c 


— table of coordinates — 


CALL JLABEL (0, 22, 
CALL JLABEL (0, 21, 
CALL JLABEL (46,18, 
CALL JLABEL (46,17, 
CALL JLABEL (46,16, 
CALL JLABEL (46,15, 
CALL JLABEL (52,14, 
CALL JLABEL (46,12, 
CALL JLABEL (46,11, 
CALL JLABEL (46,10, 
CALL JLABEL (46, 9, 


’Coordinates $’) 

’node no. X 

’For Beamlike Lattice $’) 

. $’) 

’X = coordinate along the length!’) 
’Y,Z = coordinates in the plane $’) 
’of the cross section $’) 

’For Platelike Lattice $’) 

. 

’X, Y = surface coordinates $’) 

’Z = coordinate in the $’) 


Y Z $’) 
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CALL JLABEL (52, 8, ’transverse direction $’) 

CALL JLABEL (46, 6, ’NOTE : origin must be at center $’) 

CALL JLABEL (53, 5, ’of repeating cell $’) 

DO 10 I = 1, NN 

VAL(I) = XCOORD(I) 

VAL(I+NN) = YCOORD(I) 

VAL(I+2*NN) = ZCOORD(I) 

10 CONTINUE 

HOMO = .TRUE. 

TYPEP(l) = ’R8’ 

ICOLSP(l) = 13 
FLP(l) = 10 
NROWS = NN 
NLINES = 18 

20 CALL PRTABL (0,3, NLINES, NROWS, NCOLS,ICOLSP,FLP,VAL,TYPEP, HOMO) 
CALL CKTABL (0.3, NLINES, NROWS, NCOLS.ICOLSP.FLP.VAL.TYPEP, HOMO) 
DO 30 I = 1, NN 

XCOORD(I) = VAL(I) 

YCOORD(I) = VAL(I+NN) 

ZCOORD(I) = VAL(I+2*NN) 

30 CONTINUE 
RETURN 
END 


SUBROUTINE EDSET4 

£*********************************************************************** 
C edits data set 4 - continuum strain parameters 

C********** ************************************************* ************ 

$ INCLUDE: ’ DBLE ’ 

$ INCLUDE: ’ GLOBCONS . CMN ’ 

$ INCLUDE: ’ CKCOLOR . CMN ’ 

SINCLUDE: ’ PARAM . CMN ’ 

SINCLUDE: ’DATA. CMN’ 

PARAMETER (NC0LS=3, NVAL=NCOLS*ZNSDF) 

INTEGER ICOLSP(NCOLS) , FLP(NCOLS) . VAL(NVAL) 

CHARACTER TYPEP(NC0LS)*2 
LOGICAL HOMO, ERR 
COMMON /LINES/ LINES (10) 

c — check that itruss is defined — 

IF (ITRUSS. EQ.-l) THEN 
LINES(l) = 0 

CALL ERRORS (’Please define ITRUSS in data set #1 first. Hit en 
.ter to continue $’) 

CALL WAIT 
RETURN 
END IF 

c — set defaults according to itruss — 

IF (ITRUSS. EQ.l) THEN 
NEC1 = 6 
NEC2 = 7 
NEC3 = 5 

ELSE IF ( ITRUSS. EQ. 2) THEN 
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NEC1 = 8 
NEC2 = 5 
NEC3 = 2 
ENDIF 

c — list strain parameters for itruss=l — 

CALL JLABEL (0, 22, ’Each strain parameter will be identified with 
. a number :$’) 

IF (ITRUSS.EQ. 1) THEN 

CALL JLABEL (2, 21, ’epsxo =1 gammaxyo=4 kappabr =7 epsyo=10 wbar 
=13 d(gammayzo)/dx=16$’) 

CALL JLABEL (2, 20, ’kappayo=2 gammaxzo=5 gammayzo=8 epszo=ll d(ep 
. syo)/dx=14 d(vbar)/dx =17$’) 

CALL JLABEL (2, 19, ’kappazo=3 kappat =6 psio =9 vbar =12 d(ep 
.szo)/dx=15 d(wbar)/dx =18$’) 

CALL JLABEL (55, 9, ’Defaults correspond $’) 

CALL JLABEL (55, 8, ’to Timoshenko-type $’) 

CALL JLABEL (55, 7, ’beam theory $’) 

ENDIF 

c — list strain parameters for itruss=2 — 

IF (ITRUSS.EQ. 2) THEN 

CALL JLABEL (2, 21, ’epsxo =1 kappaxo =4 gammaxzo =7 d(gam 

.maxzo)/dy=10 epszo =13$’) 

CALL JLABEL (2, 20, ’epsyo =2 kappayo =5 gammayzo =8 d(gam 

,mayzo)/dx=ll d ( epszo) /dx= 14$’) 

CALL JLABEL (2, 19, ’gammaxyo=3 2kappaxyo=6 d (gammaxzo) /dx=9 d(gam 
.mayzo)/dy=12 d(epszo)/dy=15$’ ) 

CALL JLABEL (55, 9, ’Defaults correspond $’) 

CALL JLABEL (55, 8, ’to Reissner-Mindlin $’) 

CALL JLABEL (55, 7, ’type beam theory $’) 

ENDIF 

c — edit necl ,nec2,nec3 — 

CALL JLABEL (0, 17, ’Number of strain parameters retained, NEPSC1$ 
.’) 

CALL JLABEL (0, 16, ’Number of strain parameters associated w/ neg 
. lected forces, NEPSC2 $’) 

CALL JLABEL (0, 15, ’Number of strain parameters neglected, NEPSC3 
. $’) 

ICOL = 67 


CALL PRFIELD 

(ICOL, 

17, 

2. 

NEC1 , 

’12’ , 

IBGF, 

IFGF) 

CALL PRFIELD 

(ICOL, 

16, 

2. 

NEC2, 

’12’ , 

IBGF, 

IFGF) 

CALL PRFIELD 

(ICOL, 

15, 

2, 

NEC3, 

’12’ , 

IBGF, 

IFGF) 

CALL CKFIELD 

(ICOL, 

17, 

2, 

NEC1 , 

’12’) 



CALL CKFIELD 

(ICOL, 

16. 

2. 

NEC2, 

’12’) 



CALL CKFIELD 

(ICOL. 

15, 

2, 

NEC3 , 

’12’) 




LINES (1) = NSDF 
LINES (2) = 0 

IF (NEC1.LT.1 .OR. NEC1 .GT.NSDF) THEN 

CALL ERRORS (’Please redo: hit enter and input an integer betwe 
. en 1 and $’) 

CALL WAIT 
GOTO 10 
ENDIF 
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IF (NEC2.LT. 0 .OR. NEC2.GT.NSDF .OR. 

NEC3.LT. 0 .OR. NEC3.GT.NSDF) THEN 
CALL ERRORS (’Please redo: hit enter and input an integer betwe 
.en 0 and $’) 

CALL WAIT 
GOTO 10 
ENDIF 

IF (NEC1+NEC2+NEC3 . NE . NSDF) THEN 

CALL ERRORS (’Please redo: NEPSC1+NEPSC2+NEPSC3 should equal $ 

.’) 

CALL WAIT 
GOTO 10 
ENDIF 


c — table of strain parameters — 

CALL JLABEL (0, 13, ’Choose which strain parameters are in each gr 
.oup $’) 

CALL JLABEL (0, 12, ’No. parameter retained assoc. w/negl neg 

.lected $’) 

DO 15 I * 1, NVAL 
VAL(I) = 0 
15 CONTINUE 

DO 20 I - 1, NEC1 
VAL(I) = ECU) 

20 CONTINUE 

DO 30 I = 1, NEC2 

VALU+NSDF) = ECQ+NECl) 

30 CONTINUE 

DO 40 I = 1, NEC3 

VAL(I+2*NSDF) = EC(I+NEC1+NEC2) 

40 CONTINUE 

HOMO = .FALSE. 

DO 50 I = 1, NCOLS 
TYPEP(I) = ’12’ 

ICOLSP(I) = 20+(I-l)*14 
FLP(I) = 2 
50 CONTINUE 

NROWS = NSDF 
NLINES = 10 

55 CALL PRTABL (0,2, NLINES , NROWS , NCOLS , ICOLSP . FLP , VAL , TYPEP . HOMO) 

CALL CKTABL (0,2, NLINES, NROWS, NCOLS, ICOLSP, FLP, VAL, TYPEP, HOMO) 

DO 60 1=1 ,NEC1 
EC(I) = VAL (I) 

60 CONTINUE 

DO 70 I = 1, NEC2 

ECU+NEC1) = VALU+NSDF) 

70 CONTINUE 

DO 80 I = 1, NEC3 

ECU+NEC1+NEC2) = VAL(I+2*NSDF) 

80 CONTINUE 
ERR=. FALSE. 

DO 90 I = 1, NSDF 

IF (EC(I).LT.l .OR. EC(I) .GT. NSDF) ERR= . TRUE . 

90 CONTINUE 

IF (ERR) THEN 
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LINES (1) = NSDF 
LINESC2) = 0 

CALL ERRORS (’Please redo : table numbers should be integers be 
.tween 1 and $’) 

CALL WAIT 
GOTO 55 
ENDIF 
RETURN 
END 


SUBROUTINE EDSET5 

C edits data set 5 - continuum displ. parameters 

SINCLUDE: ’ DBLE ’ 

$INCLUDE: ’ GLOBCONS . CMN ’ 

$ INCLUDE: ’ CKCOLOR . CMN ’ 

$ INCLUDE: ’ PARAM . CMN ’ 

SINCLUDE: ’DATA. CMN’ 

PARAMETER (NC0LS=2, NVAL=NCOLS*ZNDDF) 

INTEGER ICOLSP(NCOLS) , FLP(NCOLS) . VAL(NVAL) 

CHARACTER TYPEP(NCOLS) *2 
LOGICAL HOMO, ERR 
COMMON /LINES/ LINES (10) 

c — check that itruss is defined — 

IF (ITRUSS. EQ.-l) THEN 
LINES(l) = 0 

CALL ERRORS (’Please define ITRUSS in data set #1 first. Hit en 
.ter to continue $’) 

CALL WAIT 
RETURN 
ENDIF 

c — set defaults according to itruss — 

IF ( ITRUSS. EQ.l) THEN 
NUC1 = 6 
NUC2 = 6 

ELSE IF (ITRUSS. EQ. 2) THEN 
NUC1 = 5 
NUC2 = 1 
ENDIF 

c — list displacement parameters for itruss=l — 

CALL JLABEL (0, 22, ’Each displacement parameter will be identifie 
.d with a number :$’) 

IF (ITRUSS. EQ.l) THEN 


CALL 

JLABEL 

(5, 

21, 

’uo=l 

phix=4 

ubar =7 

epszo=10 

$■) 

CALL 

JLABEL 

(5, 

20, 

’ vo=2 

phiy=5 

gammayzo=8 

vbar =11 

$■) 

CALL 

JLABEL 

(5, 

19, 

’wo=3 

phiz=6 

epsyo =9 

wbar =12 

so 


CALL JLABEL (55, 9, ’Defaults correspond $’) 
CALL JLABEL (55, 8, ’to Timoshenko-type $’) 
CALL JLABEL (55, 7, ’beam theory $’) 

ENDIF 
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c — list strain parameters for itruss=2 — 

IF (ITRUSS . EQ . 2) THEN 

CALL JLABEL (5, 21. ’uo=l phix=4 $’) 

CALL JLABEL (5, 20, ’vo=2 phiy=5 $’) 

CALL JLABEL (5, 19, ’wo=3 epsz=6 $’) 

CALL JLABEL (55, 9, ’Defaults correspond $’) 

CALL JLABEL (55, 8, ’to Reissner-Mindlin $’) 

CALL JLABEL (55, 7, ’type beam theory $’) 

END IF 

c — edit nucl,nuc2 — 

CALL JLABEL (0,17, ’Number of displacement parameters retained, NUC 
.1 $’) 

CALL JLABEL (0, 16, ’Number of displacement parameters neglected, NUC 
• 2 $’) 

I COL = 50 

10 CALL PRFIELD (ICOL, 17, 2, NUC1, ’12’, IBGF, IFGF) 

CALL PRFIELD (ICOL, 16, 2, NUC2 , ’12’, IBGF, IFGF) 

CALL CKFIELD (ICOL, 17, 2, NUC1, ’12’) 

CALL CKFIELD (ICOL, 16, 2, NUC2, ’12’) 

LINES (1) = NDDF 
LINES (2) = 0 

IF (NUC1.LT.1 .OR. NUC1 .GT.NDDF) THEN 
CALL ERRORS (’Please redo: hit enter and input an integer betwe 
.en 1 and $’) 

CALL WAIT 
GOTO 10 
END IF 

IF (NUC2.LT.0 .OR. NUC2 . GT . NDDF) THEN 
CALL ERRORS (’Please redo: hit enter and input an integer betwe 
.en 0 and $’) 

CALL WAIT 
GOTO 10 
ENDIF 

IF (NUC1+NUC2 . NE . NDDF) THEN 

CALL ERRORS (’Please redo: NUC1+NUC2 should equal $’) 

CALL WAIT 
GOTO 10 
ENDIF 

c — table of displ. parameters — 

CALL JLABEL (0, 13, ’Choose which displacement parameters are in e 
.ach group $’) 

CALL JLABEL (0, 12, ’No. parameter retained neglected $’) 

DO 15 I = 1, NVAL 
VAL(I) = 0 
15 CONTINUE 

DO 20 I = 1, NUC1 
VAL(I) = UC(I) 

20 CONTINUE 

DO 30 I = 1, NUC2 

VAL(I+NDDF) = UC(I+NUC1) 

30 CONTINUE 

HOMO = .FALSE. 
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DO 50 I = 1, NCOLS 
TYPEP(I) » ’12’ 

ICOLSP(I) = 20+(I-l)*14 
FLP(I) = 2 
50 CONTINUE 

NROWS * NDDF 
NLINES = 10 

55 CALL PRTABL (0, 2, NLINES. NROWS, NCOLS, ICOLSP.FLP.VAL.TYPEP. HOMO) 
CALL CKTABL (0,2 .NLINES, NROWS .NCOLS, ICOLSP.FLP.VAL.TYPEP, HOMO) 

DO 60 I = 1, NUC1 
UC(I) = VAL(I) 

60 CONTINUE 

DO 70 I = 1, NUC2 

UCCI+NUC1) = VAL(I+NDDF) 

70 CONTINUE 
ERR=. FALSE. 

DO 90 I = 1, NDDF 

IF (UC(I).LT.l .OR. UCCI) .GT.NDDF) ERR® . TRUE . 

90 CONTINUE 

IF (ERR) THEN 
LINES (1) = NDDF 
LINES (2) = 0 

CALL ERRORS ('Please redo : table numbers should be integers be 
.tween 1 and $’) 

CALL WAIT 
GOTO 55 
END IF 
RETURN 
END 

$STORAGE : 2 

SUBROUTINE RUN 

C Runs data, does all computation 

SINCLUDE: ’ DBLE ’ 

$INCLUDE: ’ PARAM . CMN ’ 

$ INCLUDE: ’DATA. CMN’ 

C . set up . 

CCC OPEN (7 , FILE® ’ CMDAY ’ ) 

CCC WRITE (7.*) ’SETUP’ 

CALL SETUP 

C . loop over the elements . 

CALL ZEROV (KC. ZNSDF*ZNSDF) 

CALL ZEROV (PT, ZNSDF*4) 

CALL ZEROV (MC, ZNDDF*ZNDDF) 

CCC WRITE (7,*) ’LOOP OVER ELEMENTS’ 

DO 10 K - 1. NE 

C . compute element of [R] , [T] , [M] , [GAM1] , and [GAM2] . 

CCC WRITE (7,*) ’ RTMGAM , K - ’. K 
CALL RTMGAM 
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C . compute element of [KC] . [MC] . and [PT] ; add in . 

CCC WRITE (7, *) ’ KCMCPT , K = \ K 
CALL KCMCPT 
10 CONTINUE 

CCC CALL OUTMATX (KC, NSDF, NSDF, NSDF, 10, 7, ’KC : ’) 

CCC CALL OUTMATX (PT, NSDF, NSDF.NDIMT, 10, 7, ’PT : ’) 

CCC CALL OUTMATX (MC, NDDF, NDDF, NDDF, 10, 7, ’MC : ’) 

C . partition [KC] AND [MC] . 

CCC WRITE (7,*) ’PARTITION’ 

CALL PARTIT 

C . compute [C] AND [CT] . 

CCC WRITE (7,*) ’ [C] , [CT] ’ 

CALL CCT 

RETURN 

END 


SUBROUTINE SETUP 

C do all setting up : 

C 

C ... compute inverse permutations ECI and UCI for parameters 

C ... for each element, compute length 

C ... for each element, compute direction cosines X.Y.Z 

C ... for each element, compute center coordinates X,Y,Z 

$ INCLUDE: ’ DBLE ’ 

$INCLUDE : ’ GLOBCONS . CMN ’ 

SINCLUDE : ’ CKCOLOR . CMN ’ 

SINCLUDE: ’ PARAM . CMN ’ 

$INCLUDE: ’DATA. CMN’ 

PARAMETER (NC0LS=7, NVAL=NCOLS*ZNE) 

INTEGER ICOLSP(NCOLS) , FLP(NCOLS) 

REAL*8 VAL(NVAL) 

CHARACTER TYPEP(NC0LS)*2 
LOGICAL HOMO 

CALL JCLEAR 

CALL CLFIELD (0,24,80,IBGF,IFGF) 

CALL JLABEL (0,24, ’Setting up ... $’) 

C 

C . COMPUTE THE INVERSE OF (EC) AND (UC) . 

DO 2 K = 1, NSDF 
N = EC(K) 

ECI (N) = K 
2 CONTINUE 

DO 4 K = 1, NDDF 
N = UC(K) 

UCI(N) = K 
4 CONTINUE 
C 

C . LOOP OVER THE ELEMENTS . 

DO 10 K = 1, NE 



c 

c 


ccc 

c 

c 


ccc 

ccc 

c 

c 


ccc 

c 

c 


10 

c 

c 

ccc 

ccc 

ccc 

ccc 

ccc 

20 

c 

c 


30 


. FIND WHICH NODES THIS ELEMENT CONNECTS . 

I = NOP(K.l) 

J = N0P(K,2) 

WRITE (7,*) ’ELEMENT \K,’ CONNECTS NODES M.J 

. COMPUTE LENGTH . 

XX = XCOORD(J) - XCOORD(I) 

YY = YCOORD(J) - YCOORD(I) 

ZZ = ZCOORD(J) - ZCOORD(I) 

WRITE (7,*) ’ XX.YY.ZZ = ’.XX.YY.ZZ 
LEN(K) = SQRT (XX*XX + YY*YY + ZZ*ZZ) 

WRITE (7,*) ’ LEN = ’ , LEN (K) 

. COMPUTE X,Y,Z DIRECTION COSINES . 

DCOSX(K) = XX / LEN (K) 

DCOSY(K) = YY / LEN (K) 

DCOSZ(K) = ZZ / LEN (K) 

WRITE (7.*) ’ DCOSX , Y , Z= DCOSX(K) ,DCOSY(K) ,DCOSZ(K) 

. COMPUTE X,Y,Z CENTER COORDINATE . 

EX(K) = .5 * (XCOORD(I) + XCOORD(J)) 

EYCK) = .5 * (YCOORD(I) + YCOORD(J)) 

EZ(K) = .5 * (ZCOORD(I) + ZCOORD(J)) 

CONTINUE 


. WRITE OUT ELEMENT DATA . 



WRITE (7,*) ’ ’ 



WRITE (7,*) ’ ELEMENT LENGTH 

DCOSX 

DCOSY 

’ DCOSZ EX 

EY 

EZ 


WRITE (7 .20) (K.LEN(K) , DCOSX (K) .DCOSY(K) ,DCOSZ(K) .EX(K) , EY (K) , 
EZ(K) , K=1 , NE) 

FORMAT (15, E13.6, 5X, 3E13.6, 5X, 3E13.6) 


. SET UP TABLE TO PRINT . 

DO 30 I = 1, NE 

VAL(I) = LEN (I) 

VAL(I+1*NE) = DCOSX (I) 

VAL(I+2*NE) = DCOSY(I) 

VAL(I+3*NE) = DCOSZ(I) 

VAL(I+4*NE) = EX(I) 

VAL(I+5*NE) = EY (I) 

VAL(I+6*NE) = EZ(I) 

CONTINUE 
HOMO= . TRUE . 

TYPEP(l) = ’R8’ 

IC0LSP(1)= 10 
FLP(l) = 9 
NROWS = NE 
NLINES = MIN0(NE,20) 

IROW = 21-NLINES 

CALL JLABEL (19,22, ’direction cosines 
.enter $’) 

CALL JLABEL (0,21, ’Elera.no. LENGTH 1 
. Xc Yc Zc$’) 


coordinates of c 
m n 
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CALL PRTABL (0 , IROW . NLINES . NROWS , NCOLS , ICOLSP , FLP , VAL . TYPEP , HOMO) 
CALL CKTABL (0 , IROW , NLINES , NROWS , NCOLS . ICOLSP , FLP , VAL , TYPEP , HOMO) 
CALL CLFIELD (0 . 24 , 80 . IBGF , IFGF) 

CALL JLABEL (0,24,’Hit enter to continue ... $’) 

CALL WAIT 
CALL JCLEAR 

C 

C ... EXIT 
RETURN 
END 


SUBROUTINE RTMGAM 

C*** ********************************************************** ********** 

C computes for an element the following matrices : 

C 

C [R] = [L*L, M*M, N*N, L*M, L*N, M*N] 

c [10 0 ] 

C [M] = 1/6 * RHO*AREA*LEN * [2*13 13] 13= [0 1 0] 

C [ 13 2*13] . [001] 

C 

C [T] = [1, YC, ZC, YC*ZC] , BEAMLIKE 

C [T] = [1. ZC] , PLATELIKE 

C 

C [GAM1] , DEPENDS ON BEAM -OR- PLATE, CHOICE OF 
C [GAM2] , EC1 ,EC2,EC3,UC1 ,UC2 

Q*** ********************************************************** ********** 

$ INCLUDE: ’ DBLE ’ 

$ INCLUDE: ’ GLOBCONS . CMN ’ 

SINCLUDE : ’ CKCOLOR . CMN ’ 

$INCLUDE: ’ PARAM . CMN ’ 

$INCLUDE: ’DATA. CMN’ 

CALL CLFIELD (0,24, 80. IBGF, IFGF) 

CALL JLABEL (0,24, 'Computing elements number $’) 

CALL PRFIELD (28, 24, 3, K, ’12’, IBGF, IFGF) 

C 

C — COMPUTE [R] — 

R(l) = DCOSX (K) *DCOSX (K) 

R(2) = DCOSY (K) *DCOSY (K) 

R(3) = DCOSZ (K) *DCOSZ (K) 

R(4) = DCOSX (K) *DCOSY (K) 

R(5) = DCOSX (K)*DCOSZ(K) 

R(6) = DCOSY (K)*DCOSZ(K) 

CCC IF (K.EQ.l.OR.K.EQ.NE) WRITE(7,10) ’0[R], ELEMENT-’ ,K, (R(I) , 1=1 ,6) 
10 FORMAT (A, 15/ 10E13.6) 

C 

C — COMPUTE [T] — 

IF (ITRUSS.EQ. 1) THEN 
NDIMT=4 
T(l) = 1 
T(2) = EY(K) 

T(3) = EZ(K) 

T(4) = EY(K)*EZ(K) 

ELSE 
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NDIMT=2 
T(l) = 1 
T(2) = EZ(K) 

ENDIF 

CCC IF (K.EQ.l.OR.K.EQ.NE) WRITE(7,10) ’0[T] , ELEMENT 3 ' ,K, (T(I) .1-1. 
CCC . NDIMT) 

C 

C — COMPUTE [M] — 

CALL ZEROV (M, 6*6) 

IP = IPROP(K) 

VALUE = RHO (IP) *AREA(IP) *LEN (K) / 6. 

DO 20 1-1,6 

MCI, I) = VALUE*2 . 

20 CONTINUE 

M(l. 4) = VALUE 
M(2,5) = VALUE 
M(3,6) = VALUE 
M(4 , 1) = VALUE 
M(5,2) = VALUE 
M(6,3) = VALUE 

CCC IF (K.EQ.l.OR.K.EQ.NE) CALL OUTMATX (M.6,6,6, 10,7, ’ [M] : ’) 

C 

C — COMPUTE [GAM1] — 

CALL D0GAM1 

CCC IF (K.EQ.l.OR.K.EQ.NE) 

CCC .CALL OUTMATX (GAM1, 6,6, NSDF, 10,7, ' [GAM1] : ’) 

C 

C — COMPUTE [GAM2] — 

CALL D0GAM2 

CCC IF (K.EQ.l.OR.K.EQ.NE) 

CCC . CALL OUTMATX (GAM2 ,6,6, NDDF ,10,7,’ [GAM2] : ’ ) 

C 

C ... EXIT 

RETURN 
END 


SUBROUTINE D0GAM1 

c** **** **** ** * ****** ****************** ********************** ************ 
C COMPUTES ELEMENTAL MATRIX [GAM1] 

C 

C DEPENDS ON BEAMLIKE (ITRUSS-1) OR PLATELIKE (ITRUSS=2), 

C AND ON CHOICE OF EC1,EC2,EC3 

C 

C NOTE THAT EC(K) = N, WHERE THE NTH PARAMETER IS IN THE KTH POSITION 
C ECI (N)= K, 

C THEREFORE, WE NEED TO USE ECI HERE 

C********************************************************************** 

$ INCLUDE: ’ DBLE ’ 

SINCLUDE: ’GLOBCONS.CMN’ 

$ INCLUDE: ’ CKCOLOR . CMN ’ 

SINCLUDE: ’ PARAM . CMN ’ 

$ INCLUDE: ’DATA. CMN’ 

C 

C — BEAMLIKE TRUSS — 
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CALL ZEROV (GAM1, 6*NSDF) 

IF (ITRUSS.EQ. 1) THEN 
GAM1(1,ECI(D) = 1 

GAM1(1,ECI(2)) = -EY(K) 

GAM1(1 ,ECI(3)) = EZ(K) 
GAM1(1,ECI(9)) = EY(K)*EZ(K) 

GAM1(2,ECI(10)) = 1 

GAM1(2,ECI(12)) = EZ(K) 

GAM1(3.ECI(11)) = 1 

GAM1(3,ECI(13)) = EY(K) 

GAM1(4,ECI(4)) = 1 

GAM1 (4 ,ECI (6) ) = -EZ(K) 

GAM1 (4 ,ECI (7) ) = EZ(K) 

GAM1(4,ECI(14)) = EYCK) 
GAM1(4,ECI(17)) = EY(K)*EZ(K) 

GAM1(5,ECI(5)) = 1 

GAM1(5,ECI(6)) = EYCK) 

GAM1C5.ECIC7)) = EY(K) 
GAM1C5.ECIC15)) = EZ(K) 
GAM1(5,ECI(18)) = EY(K)*EZCK) 

GAM1C6.ECIC8)) = 1 

GAM1C6.ECIC12)) = EYCK) 
GAM1C6.ECIC13)) = EZ(K) 

PLATELIKE TRUSS 

ELSE 

GAMlCl.EClCl)) = 1 
GAM1C1.ECIC4)) = EZCK) 

GAM1 (2 ,ECI (2) ) = 1 

GAM1 (2 ,ECI (5) ) = EZ(K) 

GAM1(3,ECI(13)) = 1 

GAM1(4,ECI(3)) = 1 

GAM1C4.ECIC6)) = EZCK) 

GAM1(5,ECIC7)) = 1 

GAM1C5.ECIC9)) = EXCK) 

GAM1 (5 ,ECI (10) ) = EY(K) 
GAM1(5,ECI(14)) = EZ(K) 

GAM1 (6 , ECI (8) ) = 1 

GAM1(6,ECI(11)) = EXCK) 

GAM1 (6, ECI (12)) = EY(K) 
GAM1(6,ECI(15)) = EZCK) 

END IF 

. . . EXIT 
RETURN 



END 


SUBROUTINE D0GAM2 

0 ***************************************** ********** ******************** 
C COMPUTES ELEMENTAL MATRIX [GAM2] 

C 

C DEPENDS ON BEAMLIKE (ITRUSS=1) OR PLATELIKE (ITRUSS=2) , 

C AND ON CHOICE OF UC1.UC2 

C 

C NOTE THAT UC(K) = N, WHERE THE NTH PARAMETER IS IN THE KTH POSITION 
C UCI (N)= K, 

C THEREFORE, WE NEED TO USE UCI HERE 

0 *********************************** ****************************** ****** 
$ INCLUDE: ’ DBLE ’ 

$ INCLUDE: ’ PARAM . CMN ’ 

$ INCLUDE: ’DATA. CMN’ 

REAL*8 X(2) , Y(2) , Z(2) 

C 

C . FIND WHICH NODES THIS ELEMENT CONNECTS . 

I = NOP(K.l) 

J = NOP (K , 2) 

X(l) = XCOORD(I) 

X(2) = XCOORD(J) 

Y(l) = YCOORD(I) 

Y (2) = YCOORD(J) 

Z(l) = ZCOORD(I) 

Z (2) = ZCOORD(J) 

C 

C . DO FIRST USING X.YCOORD OF ITH NODE . 

C . DO SECOND USING X.YCOORD OF JTH NODE . 

CALL ZEROV (GAM2, 6*NDDF) 

DO 20 N=1 ,2 
C 

C — BEAMLIKE TRUSS — 

IF (ITRUSS.EQ.l) THEN 

GAM2(1+3*(N-1) ,UCI (1)) = 1 

GAM2(1+3*(N-1) ,UCI(5)) = Z(N) 

GAM2(l+3*(N-i) ,UCI(6)) = -Y (N) 

GAM2(1+3*(N-1).UCI(7)) = Y(N)*Z(N) 

GAM2(2+3*(N-1) ,UCI(2)) = 1 

GAM2(2+3*(N-1) ,UCI(4)) = -Z(N) 

GAM2(2+3* (N-l) ,UCI(8)) = .5*Z(N) 

GAM2(2+3*(N-1) ,UCI(9)) = Y (N) 

GAM2(2+3*(N-1) ,UCI(U))= Y(N)*Z(N) 

GAM2(3+3*(N-1) .UCI (3)) = 1 

GAM2(3+3* (N-l) ,UCI(4)) = Y(N) 

GAM2(3+3*(N-1) ,UCI(8)) = .5*Y(N) 

GAM2(3+3*(N-1) ,UCI(10))= Z(N) 

GAM2(3+3*(N-1) ,UCI(12))= Y(N)*Z(N) 

C 

C — PLATELIKE TRUSS — 

ELSE 
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GAM2(1+3*(N-1) ,UCI(1)) = 1 

GAM2 (1+3* (N-l) ,UCI(4)) = Z(N) 

GAM2(2+3*(N-1) ,UCI(2)) = 1 

GAM2(2+3*(N-1) ,UCI(5)) = Z(N) 

GAM2(3+3*(N-1) ,UCI(3)) = 1 

GAM2(3+3*(N-1) ,UCI(6)) = Z(N) 
END IF 

20 CONTINUE 
C 

C ... EXIT 
RETURN 
END 


SUBROUTINE KCMCPT 

C* ******************************************************** ************** 
C COMPUTES AN ELEMENT OF [KC] , [MC] , AND [PT] 

C 

C ELEMENT [KC] = E*AREA*LEN * ( [R] * [GAM1] )T * ([R]*[GAM1]) 

C ELEMENT [PT] = ALPHA*E*AREA*LEN * ( [R] * [GAM1] )T * [T] 

C ELEMENT [MC] = 1/OMEGA * [GAM2] T* [M] * [GAM2] 

C 

C SUM OF ELEMENTS IS DONE BY ADDING EACH ELEMENT IN TO MATRIX 
c******************** ***************************************** ********** 
$INCLUDE: ’ DBLE ’ 

$INCLUDE: ’GLOBCONS.CMN’ 

$INCLUDE: ’ CKCOLOR . CMN ’ 

$INCLUDE: ’ PARAM . CMN ’ 

$ INCLUDE: ’DATA. CMN’ 

REAL*8 RGAM1 (ZNSDF) , MGAM2(6,ZNDDF) 

REAL*8 KCTEMP(NSS) , PTTEMP(NS4) , MCTEMP(NDD) 

C 

C — ELEMENT [KC] — 

CALL MATMUL (R, GAM1 , RGAM1 , 1, 6, NSDF) 

CALL MATMUT (RGAM1 , RGAM1 , KCTEMP , 1, NSDF, NSDF) 

CCC IF (K.EQ.l.OR.K.EQ.NE) 

CCC .CALL OUTMATX (KCTEMP, NSDF, NSDF, NSDF, 10, 7, ’ELEM.KC : ’) 

IP = IPROP(K) 

VALUE = E ( IP) *AREA (IP) *LEN (K) 

CCC WRITE (7 , *) ’KC MULTIPLIER = ’, VALUE 
DO 10 I = 1, NSDF*NSDF 

KC(I) = KC(I) + VALUE* KCTEMP (I) 

10 CONTINUE 
C 

C — ELEMENT [PT] 

CALL MATMUT (RGAM1 , T, PTTEMP , 1, NSDF, NDIMT) 

CCC IF (K.EQ.l.OR.K.EQ.NE) 

CCC .CALL OUTMATX (PTTEMP, NSDF, NSDF, NDIMT, 10, 7, ’ELEM.PT : ’) 

VALUE = ALPHA (IP) *VALUE 
CCC WRITE (7 , *) ’PT MULTIPLIER = ’, VALUE 
DO 20 I = 1, NSDF*NDIMT 

PT(I) = PT(I) + VALUE* PTTEMP ( I ) 



20 CONTINUE 


C 

C — ELEMENT [MC] — 

CALL MATMUL (M, GAM2, MGAM2 , 6,6, NDDF) 

CALL MATMUT (GAM2, MGAM2, MCTEMP , 6, NDDF. NDDF) 

CCC IF (K.EQ.l.OR.K.EQ.NE) 

CCC .CALL OUTMATX (MCTEMP, NDDF. NDDF, NDDF, 10, 7. ’ELEM.MC : ’) 
DO 30 I = 1, NDDF*NDDF 

MC(I) = MC(I) + MCTEMP (I) /OMEGA 
30 CONTINUE 
C 

C ... EXIT 

100 FORMAT (A, 15/ (10E13.6)) 

RETURN 

END 


SUBROUTINE PARTIT 

0 *********************************************************************** 
C PARTITION [KC], [MC] , AND [PT] 

C 

C [KC] = [ KC11 KC12 KC13 ] [MC] = [MC11 MC12] [PT] = [PT1] 

C [ KC22 KC23 ] [MC21 MC22] [PT2] 

C [ SYM. KC33 ] [PT3] 

C 

C ONLY NEED KC11,KC12,KC21,KC22,MC11,PT1,PT2 

0 * ********************************************************************** 
$ INCLUDE: ’DBLE’ 

$ INCLUDE: ’ GLOBCONS . CMN ’ 

$ INCLUDE : ’ CKCOLOR . CMN * 

$ INCLUDE: ’PARAM.CMN ’ 

SINCLUDE: ’DATA. CMN’ 

IK(I.J) = (J-1)*NSDF + I 

IM(I , J) = (J-1)*NDDF + I 

IP(I,J) = ( J-l) *NSDF + I 

IK1(I,J) = (J-l) *NEC1 + I 
IK2(I , J) = (J-l) *NEC2 + I 
IMl(I.J) = (J-1)*NUC1 ♦ I 

CALL CLFIELD (0 , 24 , 80 , IBGF , IFGF) 

CALL JLABEL (0,24, ’Removing neglected parameters ... $’) 

C 

C — PARTITION [KC] — 

DO 10 J = 1, NEC1 
DO 10 I - 1, NEC1 

KCH(IKKI.J)) = KC(IK(I , J) ) 

10 CONTINUE 

DO 20 J = 1, NEC2 
DO 20 I - 1, NEC1 
JJ = J+NEC1 
II = I 

KC12(IK1(I,J)) = KC(IKdl.JJ)) 

KC21(IK2(J,I)) = KC(IK(JJ, II)) 

20 CONTINUE 
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DO 30 J = 1, NEC2 
DO 30 I = 1, NEC2 
JJ = J+NEC1 
II = I+NEC1 

KC22(IK2(I, J)) « KC(IK(II, JJ)) 
30 CONTINUE 


CCC 

CALL 

OUTMATX 

(KC12, 

NEC1 , 

NEC1, 

NEC2, 

10, 

7, 

’KC12 : 

’) 

CCC 

CALL 

OUTMATX 

(KC21 , 

NEC2, 

NEC2, 

NEC1 , 

10. 

7, 

’KC21 : 

’) 

CCC 

CALL 

OUTMATX 

(KCll, 

NEC1 , 

NEC1 , 

NEC1 , 

10, 

7, 

’KCll : 

’) 

CCC 

CALL 

OUTMATX 

(KC22, 

NEC2, 

NEC2, 

NEC2, 

10, 

7, 

’KC22 : 

’) 


C 

C — PARTITION [PT] — 

DO 40 J = 1. NDIMT 
DO 40 I = 1. NEC1 

PTldKld.J)) - PTdP(I.J)) 

40 CONTINUE 

DO 50 J = 1, NDIMT 
DO 50 I = 1, NEC2 
II = I+NEC1 

PT2(IK2(I, J)) = PTdPdI.J)) 

50 CONTINUE 

CCC CALL OUTMATX (PT1, NEC1, NEC1. NDIMT, 10, 7, ’PT1 : ’) 

CCC CALL OUTMATX (PT2, NEC2, NEC2, NDIMT, 10, 7, ’PT2 : ’) 

C 

C — PARTITION [MC] — 

DO 60 J = 1, NUC1 
DO 60 I - 1. NUC1 

MCll(IMld.J)) = MC(IM(I,J)) 

60 CONTINUE 

CCC CALL OUTMATX (MC11, NUC1, NUC1, NUC1, 10, 7, ’MCI 1 : ’) 

C 

C ... EXIT 

RETURN 
END 


SUBROUTINE CCT 

0 * ********************************************************************** 
C COMPUTES FINAL MATRICES 

C STIFFNESS MATRIX [CJ « 1/OMEGA * [KCll - KC12*KC22~ 1 *KC21] 

C 

C THERMAL COEFF. [CT] = 1/OMEGA * [PT1 - KC12*KC22 _1 *PT2] 

C 

C NOTES 

C ... TEMP = [KC12*KC22INV] 

C ... IF NEC2=0 , SKIP COMPUTING TEMP*KC21 AND TEMP*PT2 
C . . . USES "MATINV" TO DO MATRIX INVERSION 

C* * * * ** **** * * *********************** a|e a|e a|e afc % a|e afe 3|e ********** ****** j|e% **** ****** 

$INCLUDE: ’DBLE’ 

$ INCLUDE : ' GLOBCONS . CMN ' 

$INCLUDE: ’CKCOLOR.CMN’ 

$ INCLUDE : ’ PARAM . CMN ’ 

$ INCLUDE: ’DATA. CMN' 

COMMON /LINES/ LINES (10) 
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PARAMETER (NS2=ZNSDF*2) 

PARAMETER (MNCOLS=ZNSDF , NVAL=ZNSDF*ZNSDF) 

PARAMETER (MN2C0LS=ZNDDF , NVAL2=ZNDDF *ZNDDF ) 

REAL*8 B(ZNSDF) , IPIVOT(ZNSDF) , IWK(NS2) , TEMP(NSS) 

REAL*8 VAL(NVAL) , VAL2(NVAL) 

INTEGER ICOLSP(MNCOLS) , FLP(MNCOLS) 

INTEGER I2C0LSP(MN2C0LS) , FLP2(MN2C0LS) 

CHARACTER TYPEP(MNCOLS) *2 , TYPEP2 (MN2C0LS) *2 

LOGICAL HOMO 

IC(I.J) = (J-1)*NEC1 + I 

IM(I,J) = (J-1)*NUC1 + I 

IK2(I , J)= ( J-l) *NEC2 + I 

CALL CLFIELD (0 ,24 ,80, IBGF , IFGF) 

CALL JLABEL (0,24, ’Computing [C] and [CT] ... $’) 

C 

C . ZERO OUT MATRICES . 

CALL ZEROV (C, ZNSDF*ZNSDF) 

CALL ZEROV (CT, ZNSDF*4) 

C 

C *** ip NEC2=0 , SKIP COMPUTE KC22 INVERSE AND TEMP *** 

IF (NEC2.EQ.0) GOTO 5 
C 

C — CHECK IF ROW (& COLUMN) ARE ALL ZEROS, IF SO — 

C — REMOVE SINGULARITY WITH A ONE ON DIAGONAL 

DO 4 J = 1, NEC2 
DO 3 I = 1, NEC2 

IF (KC22 (IK2 (I , J) ) . NE . 0 . ) GOTO 4 

3 CONTINUE 

C . SINGULAR AT KC22(J,J) . 

KC22(IK2(J, J)) = 1. 

CALL CLFIELD (0,24, 80, IBGF, IFGF) 

LINES (1) = NEC1 + J 
LINES(2) = 0 

CALL ERRORS (’WARNING ! These strain parameters have no contrib 
.ution : $’) 

CALL WAIT 

CALL CLFIELD (0,24,80, IBGF, IFGF) 

4 CONTINUE 
C 

C — INVERT [KC22] — 

MROWS = 0 
I0P=0 

CALL MATINV (NEC2 , NEC2 , KC22 , MROWS , B , IOP , DETERM , ISCALE , IPIVOT , IWK) 
CCC WRITE (7 , *) ’ INVERTING KC22, DETERM = ’, DETERM 
IF (DETERM. EQ.O) STOP ’KC22 SINGULAR !’ 

CCC CALL OUTMATX (KC22, NEC2, NEC2, NEC2, 10, 7, ’KCINV : ’) 

C 

C — COMPUTE [C] — 

CALL MATMUL (KC12, KC22, TEMP, NEC1, NEC2, NEC2) 

CCC CALL OUTMATX (TEMP, NEC1, NEC1, NEC2, 10, 7, ’TEMP : ’) 

CALL MATMUL (TEMP, KC21, C, NEC1, NEC2 , NEC1) 

CCC CALL OUTMATX (C, NEC1, NEC1 , NEC1, 10, 7, ’12*22*21 : ’) 

5 DO 10 I = 1, NEC1 
DO 10 J = 1, NEC1 
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C(IC(I , J) ) = (KCIKIC(I.J)) - C(IC(I,J))) / OMEGA 
10 CONTINUE 

CCC CALL OUTMATX (C, NEC1 . NEC1, NEC1, 10, 7, ’ [C] : ’) 

C 

C — COMPUTE [CT] — 

IF (NEC2.NE.0) CALL MATMUL (TEMP, PT2, CT, NEC1, NEC2 , NDIMT) 

DO 20 I = 1, NEC1 
DO 20 J = 1, NDIMT 

CT(IC(I , J) ) = (PT1(IC(I, J)) -CTCICCI, J))) / OMEGA 
20 CONTINUE 

CCC CALL OUTMATX (CT, NEC1, NEC1, NDIMT, 10, 7, ’ [CT] : ’) 

C 

C — PRINT OUT TABLE OF [C] — 

IF (ITRUSS . EQ . 1) CALL OUTBEAM 
IF (ITRUSS. EQ. 2) CALL OUTPLAT 
DO 30 I * 1. NEC1 
DO 30 J = 1, NEC1 

VAL(IC(I.J)) = C(IC(I.J)) 

30 CONTINUE 

HOMO = .TRUE. 

TYPEP(l) = ’R8’ 

IC0LSP(1)= 4 

FLP(l) = (77- (NEC1-1) )/NECl 
FLP(l) = MIN0(FLP(1) ,15) 

NCOLS = NEC1 

NROWS = NEC1 

NLINES = NEC1 

IROW = 22-NLINES 

CALL JLABEL (0,22, ’ [C] : $’) 

CALL PRTABL (0 , IROW , NLINES , NROWS , NCOLS , ICOLSP . FLP , VAL , TYPEP . HOMO) 
CALL CLFIELD (0 , 24 , 80 , IBGF , IFGF) 

CALL JLABEL (0,24, ’Hit enter to continue ... $’) 

CALL WAIT 
CALL JCLEAR 
C 

C — PRINT OUT TABLE OF [CT] — 

DO 40 I = 1, NEC1 
DO 40 J = 1. NDIMT 

val(ic(i , j) ) = cTdcd.j)) 

40 CONTINUE 
FLP(l) = 12 
NCOLS = NDIMT 

CALL JLABEL (0,22, ’[CT] : $’) 

CALL PRTABL (0 , IROW , NLINES , NROWS , NCOLS , ICOLSP , FLP . VAL . TYPEP . HOMO) 
CALL CLFIELD (0,24,80, IBGF, IFGF) 

CALL JLABEL (0,24, ’Hit enter to continue ... $’) 

CALL WAIT 
CALL JCLEAR 
C 

C — PRINT OUT TABLE OF [MC11] — 

DO 50 I = 1, NUC1 
DO 50 J = 1, NUC1 
VAL2(IM(I , J) ) = MC11 (IM(I , J) ) 

50 CONTINUE 

HOMO = .TRUE. 


TYPEP2C1) = ’R8’ 

FLP2(1) = 12 

I2C0LSP(1)= 4 

FLP(l) = (77- (NUC1-1) )/NUCl 
FLP(l) = MIN0(FLP(1) ,15) 

NCOLS = NUC1 

NROWS = NUC1 

NLINES = NUC1 

IROW = 22-NLINES 

CALL JLABEL (0,22, ’ [m] : $’) 

CALL PRTABL (0, IROW, NLINES, NROWS, NCOLS, I2C0LSP,FLP2,VAL2,TYPEP2, 
HOMO) 

CALL CLFIELD (0,24,80, IBGF.IFGF) 

CALL JLABEL (0,24, ’Hit enter to continue ... $’) 

CALL WAIT 
CALL JCLEAR 
C 

C ... EXIT 
RETURN 
END 


SUBROUTINE OUTBEAM 

C displays output informations for beamlike lattice truss 

$ INCLUDE: ’ DBLE ’ 

$ INCLUDE: ’GLOBCONS . CMN’ 

$ INCLUDE: ’CKCOLOR.CMN’ 


CALL JLABEL (0, 22, 
.of parameters $’) 
CALL JLABEL (0, 21, 

• ~ $’) 

CALL JLABEL (0, 20, 
CALL JLABEL (0, 19, 
CALL JLABEL (0, 18, 
CALL JLABEL (0, 17, 
CALL JLABEL (0, 16, 
CALL JLABEL (0, 15. 
CALL JLABEL (0, 14, 
CALL JLABEL (0, 12, 
. of parameters $’) 
CALL JLABEL (0, 11, 

. $’) 

CALL JLABEL (0, 10, 
CALL JLABEL (4, 9, 

CALL JLABEL (4, 8. 

CALL JLABEL (4. 7, 

CALL JLABEL (4, 6, 


’ [C] for default set of strain and displ. set 


’Cll = extensional stiffness $’) 
’C22.C33 = bending stiffnesses $’) 
’C44.C55 = transverse shear stiffnesses $’) 
’C66 = torsional stiffness $’) 
’C14.C15 = extensional-shear stiffnesses $’) 
’C23 = coupling bending stiffness $’) 
’C26.C36 = bending-torsional stiffnesses $’) 


’ [CT] for default set of strain and displ. set 


’This matrix is multiplied by the vector $’) 
’( To ) $’) 

’ ( d(To)/dy ) $’) 

’ ( d(To)/dz ) $’) 

’ (d2(To)/dydz) $’) 


CALL CLFIELD (0,24,80, IBGF.IFGF) 

CALL JLABEL (0,24, ’Hit enter to continue ... $’) 
CALL WAIT 
CALL JCLEAR 
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RETURN 

END 


SUBROUTINE OUTPLAT 

C* ********************************************************************** 

C displays output information for platelike lattice truss 
£*********************************************************************** 
SINCLUDE: ’DBLE’ 

$ INCLUDE : ’ GLOBCONS . CMN ’ 

$INCLUDE: ’ CKCOLOR . CMN ’ 


CALL JLABEL (0, 22, 
.of parameters $’) 
CALL JLABEL (0, 21, 

• — - $’) 

CALL JLABEL (0. 20, 
CALL JLABEL (0, 19, 
CALL JLABEL (0. 18, 
CALL JLABEL (0, 17, 
.ffnesses $’) 

CALL JLABEL (0, 16, 
CALL JLABEL (0. 12, 
. of parameters $’) 
CALL JLABEL (0, 11, 

• - $’) 

CALL JLABEL (0, 10, 
CALL JLABEL (4. 9, 

CALL JLABEL (4, 8, 


’Cij for default set of strain and displ. set 


’ i=l,3 j=l,3 
’ i=4,6 j=4,6 
’i=7,8 j-7,8 
’i-1.3 j-4,6 


extensional stiffnesses $’) 
bending stiffnesses $’) 
transverse shear stiffnesses $’) 
bending-extensional coupling sti 


’i=4,6 j=l,3 " " " $’) 

’ [CT] for default set of strain and displ . set 


’This matrix is multiplied by the vector $’) 
’( To ) $’) 

’ (d(To)/dz) $’) 


CALL CLFIELD (0,24,80,IBGF,IFGF) 

CALL JLABEL (0,24,’ Hit enter to continue ... $’) 

CALL WAIT 

CALL JCLEAR 

RETURN 

END 


$ST0RAGE : 2 

C C 

C GENERAL MATRIX SUBROUTINES : C 

C (DOUBLE PRECISION) C 

C C 

C MATMUL ZEROV C 

C MATMUT MATINV C 

C OUTMATX C 

C- — c 


SUBROUTINE MATMUL (A,B,C,N,M,K) 

C*********************************************************************** 

MATRIX MULTIPLY: C[N,K] « A[N,M] * B[M,K] 

C*************************************************************** ******** 
$INCLUDE: ’DBLE’ 



DIMENSION A (N,M), B(M,K) , C(N,K) 
DO 10 I = 1, N 
DO 10 J = 1, K 
C(I.J) = 0. 

DO 10 L = 1, M 

C(I,J) = C(I,J) + A(I,L)*B(L, J) 
10 CONTINUE 
RETURN 
END 


SUBROUTINE MATMUT (A,B,C,M,L,N) 

C***************************************************** ****************** 

c 

C MATRIX TRANSPOSE MULTIPLY: C[L,N] = A[M,L]**T * B[M,N] 

C 

C*************************************************** ******************** 
$ INCLUDE: ’ DBLE ’ 

DIMENSION A (M.L) , B(M,N) , C(L,N) 

DO 30 K = 1, N 
DO 20 I = 1, L 
SUM = 0. 

DO 10 J = 1, M 

SUM = SUM + A(J,I)*B(J,K) 

10 CONTINUE 

C(I.K) = SUM 
20 CONTINUE 
30 CONTINUE 
RETURN 
END 


SUBROUTINE OUTMATX (K, MNROWS, NROWS. NCOLS. NFMT, IOUTF, LABEL) 
0 ******************************************************* **************** 
C OUTPUTS MATRIX [K] ON TAPE ’IOUTF’ WITH ’NFMT’ NUMBERS PER LINE 
C LABEL MUST BE 10 CHARACTERS 

Q****************** *********************************************** ****** 
SINCLUDE: ’DBLE’ 

REAL*8 K (MNROWS, NCOLS) 

CHARACTER* 10 LABEL 

WRITE (IOUTF, 10) LABEL, NROWS, NCOLS 
10 FORMAT (IX. A10/ ’ NROWS =’, 15,’ NC0LS=’,I5) 

I END = 0 

NNN = NCOLS/NFMT 

IF (MOD (NCOLS, NFMT) .NE.O) NNN = NNN+1 
DO 40 II = 1, NNN 
1ST = IEND+1 
I END = IST+NFMT-1 
IEND = MINO(IEND. NCOLS) 

WRITE (IOUTF, 50) (I, I=IST,IEND) 

DO 30 IR = 1, NROWS 

WRITE (IOUTF, 45) IR, (K(IR,J), J=IST,IEND) 

30 CONTINUE 
40 CONTINUE 
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45 FORMAT (16, 10E13.6) 
50 FORMAT (IX, 5X, 10112) 
RETURN 
END 


SUBROUTINE ZEROV (V, N) 

C*********************************************************************** 
C ZERO OUT VECTOR (V) OF LENGTH N 

Q*********************************************************************** 

$ INCLUDE: ’DBLE’ 

REAL*8 V(N) 

DO 10 I = 1, N 
V(I) = 0. 

10 CONTINUE 
RETURN 
END 


SUBROUTINE MATINV (MAX, N, A, M,B,IOP, DETERM, ISCALE, IPIVOT, IWK) 

O* *********************************** **************************************** 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


NOTE : DOUBLE PRECISION, SINGLE COMMENTED OUT WITH C* 

PURPOSE - MATINV INVERTS A REAL SQUARE MATRIX A. 

IN ADDITION THE ROUTINE SOLVES THE MATRIX 
EQUATION AX=B, WHERE B IS A MATRIX OF CONSTANT 
VECTORS. THERE IS ALSO AN OPTION TO HAVE THE 
DETERMINANT EVALUATED. 

USE - CALL MATINV (MAX, N, A .M.B.IOP, DETERM, ISCALE, IPIVOT, IWK) 

MAX - THE MAXIMUM ORDER OF A AS STATED IN THE 

DIMENSION STATEMENT OF THE CALLING PROGRAM. 

N - THE ORDER OF A, 1 . LE . N . LE . MAX . 

A - A TWO-DIMENSIONAL ARRAY OF THE COEFFICIENTS. 
ON RETURN TO THE CALLING PROGRAM, A INVERSE 
IS STORED IN A. 

A MUST BE DIMENSIONED IN THE CALLING PROGRAM 
WITH FIRST DIMENSION MAX AND SECOND DIMENSION 
AT LEAST N. 

M - THE NUMBER OF COLUMN VECTORS IN B. 

M=0 SIGNALS THAT THE SUBROUTINE IS 
USED SOLELY FOR INVERSION, HOWEVER, 

IN THE CALL STATEMENT AN ENTRY CORRE- 
SPONDING TO B MUST BE PRESENT. 

B - A TWO-DIMENSIONAL ARRAY OF THE CONSTANT 
VECTOR B. ON RETURN TO CALLING PROGRAM, 

X IS STORED IN B. B SHOULD HAVE ITS FIRST 
DIMENSION MAX AND ITS SECOND AT LEAST M. 

IOP - COMPUTE DETERMINANT OPTION. 
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I0P=0 COMPUTES THE MATRIX INVERSE AND 
DETERMINANT. 

IOP=l COMPUTES THE MATRIX INVERSE ONLY. 

DETERM- FOR I0P=0-IN CONJUNCTION WITH ISCALE 

REPRESENTS THE VALUE OF THE DETERMINANT 
OF A, DET (A) , AS FOLLOWS. 

DET (A) = (DETERM) (10**100 (ISCALE) ) 

THE COMPUTATION DET (A) SHOULD NOT BE 
ATTEMPTED IN THE USER PROGRAM SINCE IF 
THE ORDER OF A IS LARGER AND/OR THE 
MAGNITUDE OF ITS ELEMENTS ARE LARGE (SMALL) , 
THE DET (A) CALCULATION MAY CAUSE OVERFLOW 
(UNDERFLOW) . DETERM SET TO ZERO FOR 
SINGULAR MATRIX CONDITION, FOR EITHER 
I0P=1.0R 0. SHOULD BE CHECKED BY PROGRAMER 
ON RETURN TO MAIN PROGRAM. 

ISCALE - A SCALE FACTOR COMPUTED BY THE 
SUBROUTINE TO AVOID OVERFLOW OR 
UNDERFLOW IN THE COMPUTATION OF 
THE QUANTITY, DETERM. 

IPIVOT - A ONE DIMENSIONAL INTEGER ARRAY 
USED BY THE SUBPROGRAM TO STORE 
PIVOTAL INFORMATION. IT SHOULD BE 
DIMENSIONED AT LEAST N. IN GENERAL 
THE USER DOES NOT NEED TO MAKE USE 
OF THIS ARRAY. 

IWK - A TWO-DIMENSIONAL INTEGER ARRAY OF 

TEMPORARY STORAGE USED BY THE ROUTINE. 

IWK SHOULD HAVE ITS FIRST DIMENSION 
MAX. AND ITS SECOND 2. 


C REQUIRED ROUTINES- 

C 

C REFERENCE -FOX.L, AN INTRODUCTION TO NUMERICAL 

C LINEAR ALGEBRA 

C 

C STORAGE - 542 OCTAL LOCATIONS 

C 

C LANGUAGE -FORTRAN 

C 

C FORTRAN 

C LIBRARY FUNCTIONS -ABS 

C 

C LATEST REVISION - JULY 1973 -CMPB 

C 

C MATRIX INVERSION WITH ACCOMPANYING SOLUTION OF LINEAR EQUATIONS 

C******************************************* ********************************* 

$INCLUDE: ’ DBLE ’ 
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DIMENSION IPIVOT (N) , A(MAX, N) , B(MAX,N) , IWK(MAX,2) 
EQUIVALENCE (IROW.JROW), (ICOLUM, JCOLUM) . (AMAX, T, SWAP) 


c 

C INITIALIZATION 

C 

ISCALE=0 
R1 = 1.D100 
R1 = 1.E30 
R2=l . O/Rl 
DETERM=1 . 0 
DO 20 J=1 , N 
20 IPIVOT( J)=0 
DO 550 1=1, N 
C 

C SEARCH FOR PIVOT ELEMENT 

C 

AMAX=0 . 0 
DO 105 J=1,N 

IF (IPIVOT(J)-l) 60. 105, 60 
60 DO 100 K-l.N 

IF (IPIVOT(K) -1) 80. 100, 740 
80 IF (DABS ( AMAX) -DABS (A ( J , K) ) ) 85 , 100 , 100 
C* 80 IF ( ABS(AMAX) - ABS(A(J,K)))85,100,100 
85 IROW=J 
ICOLUM=K 
AMAX=A( J,K) 

100 CONTINUE 

105 CONTINUE 

IF (AMAX) 110,106,110 

106 DETERM=0 . 0 
ISCALE=0 
GO TO 740 

110 IPIVOT(ICOLUM)=IPIVOT(ICOLUM) +1 
C 

C INTERCHANGE ROWS TO PUT PIVOT ELEMENT ON DIAGONAL 
C 

IF (IROW-ICOLUM) 140, 260, 140 
140 DETERM=-DETERM 
DO 200 L=1 , N 
SWAP=A(IROW,L) 

A ( IROW , L) =A ( ICOLUM , L) 

200 A ( ICOLUM, L)=SWAP 
IF(M) 260, 260, 210 
210 DO 250 L=1 , M 
SWAP=B(IROW,L) 

B (IROW, L)=B (ICOLUM, L) 

250 B (ICOLUM, L)=SWAP 
260 IWK(I , 1)=IR0W 
IWK(I ,2)=IC0LUM 
PIVOT=A( ICOLUM, ICOLUM) 

IF(IOP.EQ.l) GO TO 321 
IF (PIVOT) 1000,106,1000 
C 

C SCALE THE DETERMINANT 
C 

1000 PIVOTI=PIVOT 

IF (DABS (DETERM) -Rl) 1030 , 1010 , 1010 
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C* IF ( ABS (DETERM) -Rl) 1030, 1010, 1010 
1010 DETERM=DETERM/R1 
I SCALE= I SCALE+ 1 

IF (DABS (DETERM) -Rl) 1060 , 1020 . 1020 
C* IF ( ABS (DETERM) -Rl) 1060, 1020, 1020 
1020 DETERM=DETERM/R1 
ISCALE=ISCALE+1 
GO TO 1060 

1030 IF (DABS (DETERM) -R2) 1040 , 1040 , 1060 
C* 1030 IF (ABS (DETERM) -R2) 1040 , 1040 , 1060 
1040 DETERM=DETERM*R1 
ISCALE=ISCALE-1 

IF (DABS (DETERM) -R2) 1050 , 1050 , 1060 
C* IF( ABS (DETERM) -R2) 1050, 1050, 1060 

1050 DETERM=DETERM*R1 
ISCALE=ISCALE-1 

1060 IF (DABS (PIVOTI) -Rl) 1090 . 1070 , 1070 
01060 IF (ABS (PIVOTI) -Rl) 1090 , 1070 , 1070 
1070 PIV0TI=PIV0TI/R1 
ISCALE=ISCALE+1 

IF (DABS (PIVOTI) -Rl) 320 , 1080 , 1080 
O IF( ABS (PIVOTI) -Rl) 320, 1080 ,1080 

1080 PIV0TI=PIV0TI/R1 
ISCALE=ISCALE+1 
GO TO 320 

1090 IF (DABS (PIVOTI) -R2) 2000 , 2000 , 320 
C*1090 IF (ABS (PIVOTI) -R2) 2000 . 2000 , 320 
2000 PIV0TI=PIV0TI*R1 
I SC ALE= I SCALE- 1 

IF (DABS (PIVOTI) -R2) 2010 , 2010 . 320 
C* IF ( ABS (PIVOTI) -R2) 2010, 2010, 320 
2010 PIV0TI=PIV0TI*R1 
ISCALE=ISCALE-1 

320 DETERM=DETERM*PIVOTI 

DIVIDE PIVOT ROW BY PIVOT ELEMENT 

321 IF (PIVOT) 330 , 106 , 330 
330 A(ICOLUM, IC0LUM)=1 .0 

DO 350 L=1 ,N 

350 A ( ICOLUM , L) =A (ICOLUM , L) /PIVOT 
IF (M) 380, 380, 360 
360 DO 370 L=1,M 

370 B (ICOLUM , L)=B( ICOLUM , L) /PIVOT 

REDUCE NON-PIVOT ROWS 

380 DO 550 Ll=l ,N 

IF(Ll-ICOLUM) 400, 550, 400 
400 T=A (LI, ICOLUM) 

A (LI, ICOLUM) =0.0 
DO 450 L=1 , N 

450 A (LI , L) =A (LI , L) -A (ICOLUM , L) *T 
IF(M) 550, 550, 460 
460 DO 500 L=1 ,M 
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500 B (LI , L) =B (LI , L) -B (ICOLUM , L) *T 
550 CONTINUE 

INTERCHANGE COLUMNS 

DO 710 1=1. N 
L=N+1-I 

IF (IWK (L , 1 ) - IWK (L , 2) ) 630 . 710 , 630 
630 JROW=IWK(L, 1) 

JCOLUM=IWK (L , 2) 

DO 705 K=1,N 
SWAP=A (K , JROW) 

A (K , JROW) =A (K . JCOLUM) 

A (K , JCOLUM) =SWAP 

705 CONTINUE 
710 CONTINUE 
740 RETURN 
END 
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Symbols 


A 

cross-sectional area of a 
member of the repeating cell 

[Cl 

matrix of stiffness coefficients 
of the simplified continuum 

Cu,C 12 , 

Cl3, • • • ,^88 

stiffness coefficients of the 
simplified continuum (see 
tables 2 and 3 and figs. 5 
and 6 ) 

{C T } 

thermal load vector of the 
simplified continuum 

di,d 2 , ■ ■ ■ ,de 

generalized displacements (see 
figs. 5 and 6 ) 

E 

elastic modulus of the material 
of a member of the repeating 
cell 

Fi,F 2 ,...,Fq 

generalized internal forces in 
the continuum beam model 
(see fig. 5) 

9 ll, 5 12 , 5 22 

partitions of the matrix [Q\ c 
(see eqs. (37)) 

\5\c 

geometric stiffness matrix of 
the continuum 

to 

geometric stiffness matrix of 
the simplified continuum 

[g (fc) ] 

geometric stiffness matrix of 
member k of the repeating cell 

[K\ 

stiffness matrix of the repeat- 
ing cell 

[K)c 

stiffness matrix of the 
continuum 

Kill Kill Ki3i 

K 22 , K 2 2 , K 33 

partitions of the matrix 
[K] c (see eqs. (14)) 

K 

kinetic energy density of the 
continuum 

L 

length of a member of the 
repeating cell 

(£,m,n)( k ) 

direction cosines of member k 

M 2 , M 3 , Mt 

bending and twisting moments 
in continuum beam model (see 
fig. 5) 

Mu, M 22 , M 12 

bending and twisting stress 
resultants in continuum plate 
model (see fig. 6 ) 


Mu, M 12 , M 22 

partitions of the matrix [M] c 
(see eqs. (15)) 

[M] c 

matrix of density parameters 
of the continuum 

[m] 

matrix of density parameters 
of the simplified continuum 
model 

mii,mi 2 , 
^13, . . . , m66 

density parameters of the 
simplified continuum (see 
tables 2 and 3 and figs. 5 
and 6 ) 

[mW] 

consistent mass matrix of 
member k of the repeating cell 

N 

axial force in beamlike lattices 
(see fig. 5) 

Ku,N 22 ,Ni 2 

extensional stress resultants 
in continuum plate model (see 
fig. 6 ) 

{ [Pt } 

thermal load vector of the 
repeating cell 

{Pt)c 

thermal load vector of the 
continuum 

{ P Tl}A P T2}i 
l P T 3 } 

partitions of the vector {Pt}c 
(see eqs. (16)) 

Ql2,Ql3 

transverse shearing forces 
in the y- and ^-directions in 
continuum beam model 

Ql,Q2 

transverse shear stress re- 
sultants in continuum plate 
model 

[J?W] 

transformation matrix whose 
entries are products of direc- 
tion cosines of member k 

Si,S 2 ,...,S$ 

stress resultants in the contin- 
uum plate model (see fig. 6 ) 

'J'O 

temperature parameter (see 
eqs. (26) and (32)) 

T 

temperature of a member of 
the repeating cell 

{T}c 

vector of temperature param- 
eters used in describing the 
continuum 

U c 

thermoelastic strain energy 
density of the continuum 
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UuU 2 

contributions to the strain 
energy of the linear and 
quadratic terms in the tem- 
perature parameters 

Uo 

isothermal strain energy 

U,V,U) 

displacement components in 
the coordinate directions 

u,v,w 

displacement parameters 
characterizing warping and 
cross-sectional distortions 

u° , v ° , w° 

displacement parameters in 
the coordinate directions 


vector of nodal displacements 
of member k of the repeating 
cell 

{u}c 

vector of displacement param- 
eters used in describing the 
continuum 

{Uclh {u c 2} 

partitions of the vector {u c } 
(see eqs. (13)) 

x,y,z 

Cartesian coordinates 

a 

coefficient of thermal expan- 
sion of a member of the re- 
peating cell 

[rl k) },[r {k) \, 

transformation matrices 
(see eqs. (6), (7), (8), and 
(34)) 

Ixyi 7z«, ~1yz 

shearing strain parameters in 
the coordinate planes 

712, 713. 723 

shearing strains in the coordi- 
nate planes 

{A} c 

vector of displacement param- 
eters and their spatial deriva- 
tives for the continuum 

{A c 1 },{A c2 } 

partitions of the vector {A} c 
(see eqs. (36)) 

£ w 

axial strain of member k of the 
repeating cell 

£ 11> £ 22> £33 

axial strains in the coordinate 
directions 



extensional strain parameters 
in the coordinate directions 

Me 

vector of strain parame- 
ters used in describing the 
continuum 

{£}(*) 

vector of strain components in 
the coordinate directions used 
in the expansion of 

{£ci } ! {^ 2 }, 
{^ 3 } 

partitions of the vector 
{e} c (see eqs. (12)) 

R 

strain parameter (see 
eqs. (25)) 

^0 K-O ^0 

5 A/ t / 1 1 

tr° 9 tr° 

’ LK, xy 

curvature changes and twist 
parameters 

p 

mass density of the material 
(see fig. 7) 

<t>xi 4>yi 

rotation components 


strain parameter (see 
eqs. (25)) 

n 

characteristic geometric 
property of the repeating 
cell of the lattice (length of 
repeating cell for beamlike 
lattices and planform area 
of repeating cell for platelike 
lattices) 

UJ 

frequency of vibration (see 
figs. 10 to 13) 

Subscripts: 


b 

bottom surface bars 

d 

diagonals 

l 

longerons 

t 

top surface bars 

V 

verticals 


Partial derivatives are denoted as follows: 

d x = d/dx, d y = d/dy, and d z = d/dz 

Superscript t denotes matrix transposition. 
Superscript ( k ) denotes member k of the repeating 
cell. 
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Table 1. Approaches for Analyzing Repetitive Lattice Structures 


Approach 

Drawbacks 

Possible remedial actions 

Direct method: 

Structure analyzed as a 
system of discrete finite 
elements 

Computationally very 

expensive for large lattices 

Use of reduction methods 
for buckling and nonlinear 
problems (refs. 43 and 44) 

Discrete field methods: 
Equilibrium and com- 
patibility equations are 
written at a typical joint 

Development and solution 

can be substantial for complex 
lattice configurations 


Taylor series expansions 
are used to develop 
differential equations 



Periodic structure approaches: 
Substructuring, combined use 
of finite elements, and 
transfer matrix methods 

Not efficient for complicated 
configurations and/or for 
transient response calculations 

Combine with substitute 
continuum approach 

Exact representation of 
stiffness of individual 
members 

Limited to beamlike lattices 

with simply supported ends or 
ring configuration with rigid 
central mast 


Substitute continuum approach 

Local deformation effects 
not accounted for 

Include local deformation 
modes in the development 


Ordinary continuum not suitable 
for lattices with rigid joints 

Use micropolar continua 
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Table 2. Stiffness Coefficients and Density Parameters for the Continuum 
Model of the Beamlike Lattice Structure Shown in Figure 7 


Coefficient 0 

Value 

Cii 

2.967 x 10 8 N 

C 14 = C 15 = C 44 = C 55 

3.126 x 10 7 N 

C22 = C33 

1.561 x 10 9 N-m 2 

C23 

9.770 x 10 7 N-m 2 

C26 = C36 = - 5^66 

-1.954 x 10 8 N-m 2 

mu = m 2 2 = m 33 

6.200 kg/m 

m44 = 2«55 = 2 m 66 

56.096 kg-m 

m 56 

1.212 kg-m 


a Cii is the extensional stiffness; C 22 and C 33 are the bending stiffnesses; C 44 and C 55 
are the transverse shear stiffnesses; Cqq is the torsional stiffness; and C 14 , C 15 , C23, C26, 
and C 36 are coupling coefficients. 


Table 3. Stiffness Coefficients and Density Parameters for the Continuum 
Model of the Platelike Lattice Structure Shown in Figure 9 


Coefficient 0 

Value 

Cn = C22 

1.682 x 10 e N/m 

Cl2 = C33 

0.4394 x 10 6 N/m 

C 14 = C 25 

1.456 x 10 6 N 

C\b = C24 = C36 

0.3802 x 10 6 N 

II 

s> 

Oi 

2.366 x 10 7 N 

C45 = Cqq 

0.6179 x 10 7 N-m 

C77 = C 8 8 

3.380 x 10 4 N/m 

mil = m22 = m 3 s 

0.1779 kg/m 2 

mi4 = m25 

0.1418 kg/m 

m44 = mss 

2.370 kg 


a C\\ and C 22 are the extensional stiffnesses in the x- and {/-directions; C 33 is the 
in-plane shear stiffness; C 44 and C 55 are the bending stiffnesses in the x- and y-directions; 
Cq 6 is the torsional stiffness; C 77 and Css are the transverse shear stiffnesses; and C 14 , C 15 , 
C 24 , C 25 , and (736 are bending-extensional coupling stiffnesses. 




Table 4. Maximum Displacements From Continuum Beam Models Compared With 
Exact Solutions for the Cantilevered Lattice Structure Shown in Figure 7 

[Number of repeating cells = 10] 


Normalized displacements 0 

Continuum model 

Exact 

Coupling terms 
neglected 

Coupling terms 
included 

u°/N, m/N 

1.685 x 10 ~ 7 

2.135 x 10 “ 7 

2.135 x 10 “ 7 

w°/Q 13 , m/N 

2.828 x 10“ 5 

3.028 x 10 ~ 5 

2.969 x 10 ~ 5 

v°/Qi 2 , m/N 

2.828 x 10~ 5 

3.028 x 10 “ 5 

2.969 x 10 “ 5 

w°/M t , N " 1 

0 

4.270 x 10 “ 7 

4.276 x 10 ~ 7 

w°/M 2 , N - 1 

8.005 x 10 - 7 

8.540 x 10 - 7 

8.442 x 10 “ 7 

v°/M 3 , N " 1 

8.005 x 10 “ 7 

8.540 x 10 ~ 7 

8.442 x 10 “ 7 


a N is the axial force; Q 13 and Q 12 are the transverse shearing forces; is the twisting moment; and 
M 2 and M 3 are the bending moments (see fig. 5). 









Table 5. Minimum Vibration Frequencies From Continuum 
Beam Models Compared With Exact Solutions for the 
Cantilevered Lattice Structure Shown in Figure 7 

[Number of repeating cells = 10] 


Mode 

Frequency from continuum model, Hz 

Exact frequency, 
Hz 

Coupling terms 
neglected 

Coupling terms 
included 

i 

3.385 

3.252 

3.280 

2 

3.385 

3.274 

3.325 

3 

13.197 

12.977 

13.136 

4 

17.059 

15.853 

16.353 

5 

17.059 

16.877 

17.172 

6 

38.794 

31.608 

31.641 

7 

38.794 

35.468 

36.116 

8 

62.185 

35.817 

37.596 

9 

62.185 

41.704 

42.996 

10 

86.209 

56.732 

59.198 


Table 6. Minimum Vibration Frequencies From Continuum 
Plate Models Compared With Exact Solutions for the 
Cantilevered Double-Layered Grid Shown in Figure 9 


Mode 

8 L x 8 L grid 

16L x 16L grid 

Continuum 

frequency, 

Hz 

Exact 

frequency, 

Hz 

Continuum 

frequency, 

Hz 

Exact 

frequency, 

Hz 

1 

1.253 

1.248 

0.3887 

0.3878 

2 

2.064 

2.087 

0.7152 

0.7178 

3 

4.260 

4.249 

1.703 

1.700 

4 

4.768 

4.792 

2.182 

2.189 

5 

5.527 

5.573 

2.346 

2.345 

6 

6.303 

6.298 

2.387 

2.391 

7 

8.273 

8.347 

3.436 

3.455 

8 

8.377 

8.396 

3.670 

3.664 

9 

9.090 

9.144 

3.979 

3.984 

10 

9.856 

9.809 

4.360 

4.354 
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Double-layered tetrahedral platelike lattice Continuum plate model 

Figure 1. Continuum models for beamlike and platelike lattice structures. 








Figure 2. Local deformations in planar lattice trusses subjected to axial loading P. 



Truss member (pin joints) 

Beam member (rigid joints) 

Deformation 

pattern 





i L j 

I - * L ** 

" '• 

W. 

1 

<. 1 ■' 1 j 

w. 

!. 

J 



Rotations 

V = (w. - w. ) 

V = — (\ w. - w. ) 

L J i 

0j,0j (joint rotations) 

Appropriate 

continuum 

Ordinary 

(displacement field only) 

Micropolar 

(independent displacement 
and rotation fields) 


Figure 3. Deformation patterns for pin-jointed and rigid-jointed one-dimensional members. 
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Phase 1 
3-D 

Continuum 


Phase 2 
1- or 2-D 
Continuum 


Phase 3 
Simplified 
1- or 2-D 
Continuum 



Figure 4. Steps involved in the generation of the stiffness coefficients for simplified one- and two-dimensional 
continuum models of beamlike and platelike lattice structures. 
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and associated stiffness and density parameters. 
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A t = 80 x 10" 6 m 2 , E = 71.7 x 10 9 Pa 

A b = 50x10' 6 m 2 

Ajj = 10 x 10 ' 6 m 2 

A v = lOx 10‘ 6 m 2 . p = p b =p d = p v 

L = 7.5 m = 2768 kg/m 3 

h = 7.5 m 



Cross 

section 

area 

Length 

Desig. 

Top surf, 
bars 

A t 

L,f2L 


Bottom 
surf, bars 

A b 

l,/2T 


Verticals 

A v 

h 


Diagonals 

A d 

d 



Figure 9. Double-layered hexahedral platelike lattice truss. 
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Mode number 
(a) Ten-bay truss. 



Mode number 
(b) Thirty-bay truss. 

Figure 10. Accuracy of minimum frequencies from continuum beam models for the cantilevered beamlike lattice 
structure shown in figure 8. 
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co 3 = 1 2.98 (1 3.14) co 4 = 1 5.85 (1 6.35) 



co 5 = 16.88(17.17) to 6 = 31 .61 (31.64) 



® 7 = 35.47 (36.12) co 8 = 35.82 (37.60) 



w 9 =41.70(43.00) 



<o 10 = 56.73 (59.20) 


Figure 11. Vibration mode shapes for the cantilevered lattice 10-bay truss shown in figure 8 
parentheses are exact frequencies (in hertz); other numbers are continuum beam frequencies. 


Numbers in 
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Mode number 

(a) 8 L x 8 L grid. 


1.05 



Mode number 
(b) 16L x 16L grid. 


Figure 12. Accuracy of minimum frequencies from continuum plate models for the cantilevered platelike lattic 
truss shown in figure 10. 
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C0 1 = 1 .253 (1 .248) 







Figure 13. Vibration mode shapes for the cantilevered hexahedral grid shown in figure 10. Numbers 
parentheses are exact frequencies (in hertz); other numbers are continuum plate frequencies. 
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